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data.  In  the  long  term,  we  will  design  custom  arrays  for  use  in  clinical  practice.  Genes  will  be  further  studied  using  cellular  and 
molecular  methods,  and  their  role  as  therapeutic  targets  explored. 
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Endocrine  therapy  is  often  the  least  toxic  and  most  effective  treatment  for  hormone  receptor  positive  invasive 
breast  cancer.  Such  therapy  includes  antiestrogens  (tamoxifen,  fulvestrant)  and  aromatase  inhibitors 
(anastrozole,  letrozole,  exemestane).  Tamoxifen  (TAM)  increases  disease  free  and  overall  survival  in  the 
adjuvant  setting,  reduces  the  incidence  of  estrogen  receptor  positive  disease  (ER+;  unless  otherwise  noted 
ER=ERa)  in  high-risk  women,  and  reduces  the  rate  of  bone  loss  secondary  to  osteoporosis  in  postmenopausal 
women  [1,2].  Aromatase  inhibitors  are  effective  only  in  the  absence  of  functioning  ovaries  -  TAM  can  be  used 
regardless  of  menopausal  status.  Recent  studies  suggest  that  anastrozole  may  be  superior  to  TAM  in  the 
adjuvant  treatment  of  postmenopausal  women  with  ER+  breast  cancer;  other  studies  report  higher  overall 
response  rates  with  letrozole  (LET)  vs.  TAM  as  first  line  therapy  in  the  metastatic  setting.  Thus,  a  recent 
controversy  in  the  management  of  patients  with  ER+  disease  is  whether  an  aromatase  inhibitor  or  TAM  should 
be  given  as  first  line  endocrine  therapy  [3-9]. 

In  this  new  Clinical  Translational  Research  award,  we  will  build  classifiers  that  accurately  separate 
antiestrogen  sensitive  from  antiestrogen  resistant  breast  tumors  and  begin  to  assist  in  the  direction  of  specific 
endocrine  treatments  (antiestrogen  vs.  aromatase  inhibitor)  to  individual  patients.  We  hypothesize  that 
endocrine  responsiveness  is  affected  by  a  gene  network,  rather  than  the  activity  of  only  one  or  two  genes  or 
signaling  pathways  [10-12].  Since  the  key  components  of  such  a  network  are  unknown,  we  must  study  10,000s 
of  genes.  We  will  use  Affymetrix  GeneChips.  We  will  not  identify  mutational  events,  the  presence  of  mRNA 
splice  variants,  or  post-translational  protein  modifications.  However,  these  factors  have  major  effects  on  the 
transcriptome  and  their  "footprints"  should  be  identified  by  expression  microarrays. 


Body 

Overview:  We  will  build  classifiers  that  separate  antiestrogen  sensitive  from  antiestrogen  resistant  breast 
tumors  and  begin  to  assist  in  the  direction  of  specific  endocrine  treatments  (antiestrogen  vs.  aromatase  inhibitor) 
to  individual  patients.  To  achieve  this  goal,  and  consistent  with  a  CTR  award,  we  will  complete  a  4-year, 
prospective,  neoadjuvant  study  with  Letrozole  (LET)  or  TAM  as  the  only  systemic  therapy.  We  will  obtain 
molecular  profiles  from  Affymetrix  GeneChips  and  further  develop  and  apply  our  innovative  bioinformatic  and 
biostatistic  methods  to  explore  these  high  dimensional  data  sets  and  build/validate  new  classifiers.  A  more 
accurate  predictor  of  endocrine  responsiveness  would  have  widespread  clinical  use,  allowing  women  and 
physicians  to  make  more  individualized  and  appropriate  treatment  decisions.  For  example,  patients  with  tumors 
predicted  to  be  resistant  to  antiestrogens  and/or  aromatase  inhibitors  would  be  strong  candidates  for  an  early 
intervention  with  cytotoxic  chemotherapy. 

In  most  predictive/prognostic  marker  studies  investigators  focus  on  a  single  factor  and  whether  they 
obtain  a  p-value  that  reaches  conventional  statistical  significance.  Our  approach  is  different  because  we  will 
determine  whether  we  can  find  joint  gene  subsets  that  can  separate  patients  into  sufficiently  distinct  groups  that 
should  differ  in  their  treatment.  We  will  (1)  analyze  >33,000  genes  on  retrospective  and  prospective  material, 
(2)  apply  new  biostatistical  and  bioinformatic  methods  to  identify  ~40  potentially  informative  "biomarkers,"  (3) 
build  neural  network  and  biostatistical  model  classifiers,  (4)  evaluate  the  joint  discriminant  power  of  selected 
genes  concurrently  rather  than  as  single  biomarkers,  (5)  focus  on  prediction  for  individual  patients  where  the 
assessment  of  a  p-value  is  less  important  than  the  classification  rate  of  our  predictors,  (6)  validate  the  classifiers 
in  independent  data  sets,  and  (7)  explore  the  ability  of  predictors  to  refine  the  targeting  of  specific  endocrine 
therapies. 

Evidence  has  begun  to  accumulate  suggesting  that  an  aromatase  inhibitor  might  be  a  more  effective  first 
line  endocrine  therapy  for  some  breast  cancer  patients  than  the  current  standard  of  care  (Tamoxifen).  These  data 
have  generated  considerable  interest  and  controversy,  in  part  because  unlike  TAM,  there  are  no  long  term 
studies  with  aromatase  inhibitors  where  definitive  survival  data  are  available.  Our  study  could  provide  new  and 
innovative  insights  into  how  to  approach  the  more  effective  targeting  of  specific  endocrine  therapies  to 
individual  patients. 
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Specific  Aims  (from  the  original  application) 

We  will  complete  two  clinical  studies  and  collect  gene  expression  profiles  from  which  to  build  predictors  of 
endocrine  responsiveness.  Predictors  will  be  built  in  Specific  Aim  2  and  validated  in  Specific  Aim  3. 

Aim  1:  Clinical  Studies  -  Clinical  Studv-1  (retrospective)  is  of  pretreatment,  single,  frozen  samples  where  we 
will  compare  the  molecular  profiles  of  tumors  that  recurred  on  TAM  with  those  of  tumors  that  did  not  recur. 
Each  resistant  sample  is  matched  with  a  TAM  sensitive  sample  by  age,  stage,  and  duration  of  follow-up.  We 
also  have  further,  single  (unmatched),  frozen  samples  from  patients  already  progressing  on  TAM.  Clinical 
Studv-2  is  a  prospective  study  of  breast  tumor  samples  from  patients  treated  with  neoadjuvant  TAM  or  LET. 

Aim  2:  We  will  develop  and  apply  novel  bioinformatics  and  biostatistics  to  discover  gene  subsets  that  define  the 
molecular  differences  between  endocrine  sensitive  and  resistant  breast  tumors.  These  genes  will  be  used,  in 
combination  with  established  predictive/prognostic  factors,  e.g.,  ER,  PgR,  stage,  to  build  innovative  classifiers 
that  can  better  predict  an  individual  tumor’s  endocrine  responsiveness. 

Aim  3:  We  will  test,  optimize,  and  validate  the  performance  of  the  classifiers  from  Aim  2  in  retrospective 
studies  of  human  breast  tumors.  We  will  measure  each  gene  individually  by  IHC,  in  situ  RNA  hybridization 
(ISH),  or  real  time  PCR  (RT-PCR). 


Key  Research  Accomplishments 

Progress  on  the  clinical  goals  for  Year-1  was  greatly  delayed  because  of  the  time  taken  to  obtain  DOD  approval 
of  our  preexisting  institutionally  approved  IRBs  at  Georgetown  University  and  at  the  University  of  Edinburgh. 
All  institutionally  approved  protocols  and  requested  material  were  submitted  to  the  DOD  in  July  2004; 
additional  information,  requested  by  the  DOD  several  months  later,  was  submitted  in  November  2004.  We  did 
not  receive  final  approval  to  proceed  with  the  clinical  studies  until  March  2005.  Much  of  this  delay  seems  to 
have  been  unavoidable  and  due,  in  part,  to  personnel  changes  at  the  DOD  (within  USAMRMC).  While  this  has 
likely  left  us  behind  schedule  in  recruitment  to  the  prospective  studies,  we  were  able  to  proceed  with  the 
informatics  (algorithm  development  and  optimization)  and  infrastructure  development  (database  development 
and  installation).  We  also  used  this  time  to  obtain  retrospective  material  from  the  University  of  Edinburgh  (we 
obtained,  annotated,  and  stored  this  material  but  could  not  use  the  samples  until  we  had  formal  final  approval 
from  the  DOD),  and  to  further  optimize  our  microarray  procedures  and  complete  modifications  to  our  standard 
operating  procedure  documents.  We  also  completed  some  of  the  studies  presented  in  our  preliminary  data  and 
submitted  these  for  publication.  The  in  silico  tissue  heterogeneity  correction  method  described  in  the  application 
was  developed  sufficiently  and  submitted  for  publication  -  this  was  sent  out  for  review  by  Nature  Methods  but 
we  do  not  know  the  final  outcome.  Publications  supported  since  this  award  are  listed  under  “Reportable 
Outcomes”  (below)  and  constitute  some  of  our  major  accomplishments  in  the  first  year. 

Other  key  research  accomplishments  are  related  to  the  goals  set  towards  the  establishment  of  our  Center 
application;  these  are  more  closely  related  to  the  Statement  of  Work  and  are  presented  in  this  context  below. 

Statement  of  Work  (from  the  original  application) 

•  TASK  1.  Array  breast  tumor  samples  from  Clinical  Studies  1  (retrospective)  and  2  (prospective):  (1-48 
Months) 

To  perform  this  task  we  will  obtain  breast  tumor  samples  and  clinical  information  from  University  of 
Edinburgh,  collect  and  quality  test  RNA  using  validated  tissue  acquisition  and  processing  protocol,  and  array 
RNA  samples  on  oligonucleotide  chips  (i.e.,  U133A  Affymetrix  GeneChips).  Please  note  that  we  originally 
described  analyses  of  approximately  12,000  genes  in  each  sample  and  now  indicate  that  we  will  measure 
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almost  3-times  as  many  genes.  The  increase  is  possible  because  Affymetrix  improved  their  technology  and  now 
produce  single  chips  with  40,000  probe  sets  representing  39,000  transcripts,  of  which  33,000  are  well- 
substantiated  genes.  The  cost  of  these  chips,  which  essentially  represent  the  probe  sets  previously  included  on 
two  chips  (U133A  and  U133B),  is  the  same  as  the  original  U133A  GeneChips  described  when  the  application 
was  submitted.  Since  we  were  unable  to  start  arraying  in  year  1  (see  below),  this  proved  very  fortunate,  since  it 
greatly  increases  the  power  of  our  study  to  detect  meaningful  predictive  patterns  and  genes  or  networks 
associated  with  the  clinical  outcomes. 

We  have  recently  begun  extracting  and  assessing  RNA  from  the  samples  we  have  received  to  date.  If  the 
high  quality  of  the  initial  specimens  is  apparent  in  subsequent  material,  we  should  be  able  to  start  arraying  the 
retrospectively  collected  material  within  the  next  few  months.  Thus,  we  would  hope  to  remain  on  track  in  Year 
2  and  begin  generating  our  first  series  of  data  for  analysis. 

Note:  minor  change  to  SOW  -we  will  use  the  Affymetrix  U133  Plus  2.0  GeneChips  (not  U133A). 


•  TASK  2.  Store,  process,  and  train/optimize  classifiers  from  gene  expression  microarray  data:  (1-48 
Months) 

To  perform  this  task  we  will  modify  MIAME  Compliant  DSAD  Database  that  can  be  readily  exported  to 
National  Cancer  Institute  Center  for  Bioinformatics.  We  will  submit  collected  de-identified  clinical  information 
into  the  DSAD  database  and  process  gene  expression  data  with  “in  house”  state-of-the-art  algorithms  (we  will 
also  further  develop  and  optimize  these  algorithms  throughout  this  award  period).  For  the  initial  studies  we  will 
train/optimize  initial  neural  network  RNA  classifier  (MLP),  the  final  classifier  for  the  microarray  data  will  be 
built  when  we  have  completed  arraying  all  samples. 

We  have  made  significant  progress  on  addressing  this  task  and  are  likely  ahead  of  schedule,  largely  as  a 
consequence  of  our  involvement  in  the  National  Cancer  Institute  Center  for  Bioinformatics  (NCICB)  caBIG 
project.  The  PI  (Dr.  Clarke)  leads  the  Lombardi  Comprehensive  Cancer  Center’s  caBIG  team  and  we  have  been 
actively  involved  in  the  development  of  caArray  (NCICB’s  grid-enabled,  MIAME  compliant,  microarray 
database).  We  chose  to  adopt  this  database,  rather  than  continue  working  on  our  own  DSAD,  because  caArray  is 
further  developed  and  already  has  most  of  the  functions  we  require.  We  now  have  a  working  installation  of 
caArray  and  are  already  able  to  store  and  exchange  microarray  data  in  our  testing  of  this  system.  The  caBIG 
program  is  open  source-open  access  and  is  widely  supported  by  NCICB  and  teams  of  collaborating  scientists  at 
other  Cancer  Centers  across  the  country.  We  also  have  found  the  NCICB  team  highly  responsive  when  we 
identify  bugs  or  problems  with  the  software.  We  anticipate  that  continued  collaboration  through  the  caBIG 
community  will  prove  a  more  cost  and  time  efficient  approach  to  developing  some  components  of  the  research 
infrastructure  described  in  the  original  application.  It  is  our  intent  to  build  any  additional  components  in  a 
manner  consistent  with  the  guidelines  established  by  the  caBIG  community,  since  this  will  likely  ensure  long 
term  viability  and  compatibility  of  our  database  structure. 

Note:  minor  change  to  SOW  -  we  will  now  use  the  caArray  database  rather  than  the  DSAD  as  we  initially 
proposed. 


•  TASK  3.  Retrain/reoptimize  classifiers  using  IHC  data  from  Series  1  (Archival  Tissues)  and  Series  2 
(Scottish  Adjuvant  TAM  Trial)  for  Validation:  (Months  24-48): 

To  perform  this  task  we  will  obtain  clinical  information  and  breast  tumor  samples  from  the  University  of 
Edinburgh  (formalin  fixed/paraffin  embedded).  We  will  rank  and  prioritize  selected  joint  genes  from  RNA 
classifiers  built  and  optimized  in  TASK  2  (above)  and  retrain/reoptimize  the  initial  neural  network  IHC 
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classifier  (MLP).  Finally,  we  will  validate  IHC  classifier  on  independent  data  sets  (data  sets  not  used  to  build 
and  train  the  MLP  classifiers). 

We  will  not  be  able  to  start  this  task  on  the  timeframe  as  initially  proposed  here  because  of  the  delays  in  getting 
approval  to  work  with  the  clinical  specimens.  However,  we  hope  to  be  able  to  recover  time  in  the  coming  12 
months  and  still  finish  these  tasks  within  the  period  of  funding.  We  have  completed  the  preliminary  data  on  the 
use  of  IHC  and  tissue  microarrays  to  explore  an  initial  gene  set  implicated  in  antiestrogen  resistance.  We  also 
have  performed  and  published  initial  mechanistic  studies  on  some  of  these  genes,  consistent  with  the  goals  of 
Task  3,  since  we  anticipate  obtaining  additional  data  further  implicating  these  genes  in  Task  1. 


Reportable  Outcomes 
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*We  include  in  the  appendix  reprints  of  those  papers  that  are  already  published.  Manuscripts  cited  as  “in  press” 
will  be  included  in  the  next  annual  report,  once  reprints  are  available.  We  do  not  include  the  abstracts,  but  can 
do  so  if  requested  by  either  the  reviewers  and/or  DOD.  Several  other  manuscripts  related  to  our  bioinformatic 
methods  also  are  submitted  and  in  preparation;  these  will  be  included  as  appropriate  in  later  annual  reports 

Conclusions 

We  have  made  good  progress  on  the  research  infrastructure  goals  and  in  the  development  or  optimization  of  the 
methods  needed  for  data  analysis.  We  also  have  completed  and  published  most  of  the  data  presented  as 
preliminary  data  in  the  initial  application.  The  clinical  studies  were  held  up  by  an  unexpectedly  long  delay  in 
obtaining  final  approval  for  our  existing  protocols  but  this  is  now  taken  care  of  and  we  are  poised  to  begin 
analysis  of  our  first  series  of  breast  cancer  specimens.  We  are  confident  that  we  can  maintain  a  high  level  of 
productivity  and  remain  on  target  to  complete  the  studies  by  the  end  of  the  funding  period. 
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Abstract.  Multilayer  perceptrons  offer  an  integrated  procedure  for  feature 
extraction  and  Bayes  classification  by  learning  the  decision  boundary.  Its 
feed-forward  autoassoclative  architecture  can  also  be  used  to  construct 
subspaces  in  a  supervised  or  unsupervised  model  [1].  On  the  other  hand, 
multiclass  linear  discriminant  analysis  provides  a  multivariate  prediction  by 
estimating  the  density  function.  Its  linear  subspaces  obtained  by  the  weighted 
Fisher  criteria  under  a  standard  finite  normal  mixture  model  retain  most 
closely  the  intrinsic  Bayes  separability  [2J.  Here  we  show  a  twofold  connection 
between  multilayer  perceptrons  and  linear  discriminant  analysis  using 
discriminatory  component  analysis.  This  theoretical  observation  immediately 
suggests  a  possible  clustering-model  supported  optimization  mechanism  for 
multilayer  perceptrons:  the  weights  between  the  input  and  hidden  layers  are 
related  to  eigenvectors  of  the  weighted  Fisher  scatter  matrix,  the  number  of 
the  hidden  layer  neurons  is  justified  by  the  corresponding  significant 
eigenvalues,  and  the  weights  connected  to  the  output  neurons  are  obtained 
from  the  centers  of  the  classes  in  the  extracted  feature  subspaces. 


INTRODUCTION 

Multilayer  perceptrons  (MLP)  offer  an  integrated  procedure  for  feature  extraction 
and  Bayes  classification  by  learning  the  decision  boundary.  Its  feed-forward 
autoassociative  architecture  can  also  be  used  to  construct  nonlinear  subspaces  in  a 
supervised  or  unsupervised  model  [1,3].  The  output  of  each  hidden  layer  may  be 
interpreted  as  a  set  of  new  features  presented  to  the  output  layer  for  classification 
[1].  Such  feature  extraction  actually  has  various  functions.  For  example,  one 
typical  utility  is  to  perform  dimensionality  reduction  [1],  Another  popular 
application  is  to  transform  a  nonlinearly  separable  problem  into  a  linear 
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classification  [3].  When  both  functions  are  required,  multiple  hidden  layers  may  be 
needed  [1,3]. 

On  the  other  hand,  multiclass  linear  discriminant  analysis  provides  a 
multivariate  prediction  by  estimating  the  density  function.  Its  subspaces  obtained 
by  the  weighted  Fisher  criteria  (wFC)  under  a  standard  finite  normal  mixture 
(SFNM)  model  retain  most  closely  the  intrinsic  Bayes  separability  [2].  It  can  be 
shown  that  the  determination  of  the  linear  dimension  reduction  (LDR)  transform  is 
equivalent  to  finding  the  maximum-likelihood  parameter  estimates  of  a  SFNM 
model  [2,4].  This  motivates  an  exploration  of  the  connections  between  MLP  and 
LDR.  A  natural  hypothesis  is  that  the  class  labels  are  used  as  targets  during 
supervised  training,  forcing  the  outputs  of  the  hidden  layer  to  capture  the  most 
discriminatory  components  or  subspaces,  so  as  to  maximize  the  classification 
performance  of  the  output  layer  [1].  Furthermore,  the  output  neuron  in  MLP  plays 
a  similar  role  as  the  perceptron.  It  bears  a  close  resemblance  to  Fisher’s  linear 
discriminant  analysis  (LDA)  [3],  and  has  been  generalized  to  multiple  classes  [2,5], 

Based  on  these  theoretical  observations,  we  therefore  wish  to  suggest  an 
optimization  mechanism  for  multilayer  perceptrons.  The  issues  at  hand  include 
initialization  of  the  weights  between  input,  hidden,  and  output  layers;  and  the 
determination  of  the  numbers  of  hidden  and  output  neurons.  We  demonstrate  the 
principle  of  the  method  on  both  simulated  and  real-world  data  sets. 


THEORY  AND  METHOD 

When  there  are  more  than  two  classes,  the  task  is  intrinsically  a  nonlinear 
classification  problem.  Our  initial  discussion  will  focus  on  multiclass  cases  where 
all  class-pairs  are  linearly  separable,  and  we  will  then  propose  our  strategies  to 
more  general  nonlinearly  separable  cases. 

Linear  Dimension  Reduction  and  Feature  Extraction 

Consider  a  given  mo-dimensional  input  t-space,  multiclass  LDR  is  to  search  for  a 
linear  transformation  W  that  reduces  the  dimensionality  to  a  lower  mt- 
dimensional  feature  x-space  (mi«mo)>  while  preserving  a  maximum  amount  of 
discrimination  information.  Since  it  is  too  complex  to  use  the  Bayes  error  directly 
as  a  criterion,  the  most  well-known  technique  is  LDA  that  determines  the 
projection  matrix  such  that  the  Fisher  criterion  of  total  scatter  versus  average 
within-class  scatter  is  maximized  [1].  It  has  been  shown  that  when  there  are  more 
than  two  classes  (Ko>2),  the  Fisher  criterion  that  uses  the  squared  Mahanalobis 
distance  between  the  classes  in  the  dimension  reduced  subspace  is  suboptimal  with 
respect  to  classification,  since  large  between-class  distances  are  overemphasized 
and  the  resulting  projection  preserves  the  distances  of  already  well-separated 
classes,  causing  a  large  overlap  of  neighboring  classes  [2]. 

By  a  decomposition  of  individual  class  pairs,  the  wFC  was  recently  proposed 
in  order  to  improve  upon  LDA  in  which  the  weighting  for  the  contribution  is 
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derived  from  an  attempt  to  approximate  the  Bayes  error  for  pairs  of  classes,  taking 
a  sum  of  the  following  general  form  [2] 

^rc(W)  =  fX/rMAJtrace^SiXw)  (1) 

/.*4l 


where  ir*  and  are  the  prior  probabilities  of  classes  k  and  /,  S,„.=2jrakC,*  is  the 
pooled  within-class  scatter  matrix,  S^r(Pti'lb[)(mrPit)T  is  the  between-class 
scatter  matrix,  and  a>(A*/)  is  the  weighting  function  that  depends  on  the 
Mahalanobis  distance  Ajpftptt-pt/fSh^ptrpu)]1'2  between  the  classes  k  and  l, 
with  class  mean  vector  p,  and  covariance  matrix  C(.  Furthermore,  the  weighting 
function  that  approximates  the  Bayes  error  rate  between  the  classes  via  criterion 
(1)  can  be  shown  to  take  an  expression  of  the  following  form  [2] 


where  A*r=||pt*-pt/||  comes  down  to  the  ordinary  Euclidean  distance  between  class 
means,  and  erf(.)  is  the  standard  Gaussian  error  function  [2]. 

Finding  a  solution  W  that  maximizes  such  a  criterion  (1)  comes  down  to 
determining  an  eigenvalue  decomposition  of  the  matrix 

S'<L  £  JSw  (3) 


and  constructing  W  by  taking  the  /»]  eigenvectors  corresponding  to  the  m,  largest 
eigenvalues  [1],  namely  discriminatory  component  analysis  (DCA).  Although  both 
St*,'1  and  St*/  are  symmetric,  however,  matrix  (3)  is  generally  not  symmetric, 
resulting  in  possibly  complex-valued  eigenvalues  and  eigenvectors.  An  elegant  yet 
simple  alternative  is  to  determine  an  eigenvalue  decomposition  of  the  generalized 
between-class  scatter  matrix 


*=l  /=,**! 


and  constructing  the  interim  LDR  matrix  W’  similarly  as  aforementioned.  Notice 
that  matrix  (4)  is  always  symmetric.  It  can  be  shown  that  W-WS^'1'2  maximizes 
the  original  criterion  (1). 

Having  determined  LDR  matrix  W,  the  dimension  reduced  feature  subspaces 
become  x=WT(t-b,o),  where  b(0  is  the  global  center  of  the  data  set,  and  x/=w/(t- 
bw)  for  y=l,  ...,  in |.  On  the  other  hand,  the  output  of  the  hidden  layer  in  MLP  is 
computed  as  xr=<p(\\jy t-bj)  where  ip{.)  is,  often  nonlinear,  the  activation  function. 
The  similarity  between  LDR  and  feature  extraction  immediately  suggests  that  the 
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column  vectors  of  the  LDR  matrix  W  can  be  used  to  initialize  the  weights  between 
the  input  and  hidden  layer  of  an  MLP  with  a 
bias  6/=w/b,o,  except  for  a  nonlinear  scaling 
of  each  new  feature.  Accordingly,  the 
number  of  the  hidden  neurons  m  i  shall  be 
justified  by  the  number  of  the  significant 
eigenvalues  of  matrix  (4).  It  has  been 
theoretically  shown  that  minimization  of  the 
Bayes  error  with  respect  to  the  synaptic 
weights  of  the  MLP  is  equivalent  to 
maximizing  the  criterion  (1)  that  is 
determined  entirely  by  the  hidden  neurons  of 
the  MLP  [3], 


Linear  Discriminant 
Multiclass  perceptrons 


Analysis  and 


Figure  1.  Decision  making  by  a 
perceptron  of  2-output  neurons. 


As  primarily  a  two-class  classifier,  LDA  describes  a  linear  transformation  from  an 
mi-dimensional  problem  to  a  one-dimensional  problem  [3],  Consider  a  variable 
y=v/r*-b,  the  LDA  is  defined  by 


JO): 


:  ™TSlt;W 
wTSI„w 


(5) 


that  is  known  as  the  generalized  Rayleigh 
quotient.  The  solution  that  maximizes 
J(w)  is  simply  w=S™ which  is 
also  a  generalized  eigenvalue  problem  and 
referred  to  as  Fisher’s  linear  discriminant 
[3]. 

As  aforementioned,  the  output  neuron 
in  the  MLP  functions  similarly  as  the 
perceptron  that  bears  a  close  resemblance 
to  the  Bayes  classifier,  and  has  been 
generalized  to  multiple  classes  [3]. 
Specifically,  the  output  neuron  in  the 
MLP  is  computed  as  y,=$)(w,Tx-bi),  for 
i=  1,  ....  m2.  Consider  a  two-class  case 


Figure  2.  Decision  making  by 
perceptron  of  3-output  neurons. 


with  a  linear  activation  function,  we  have  y=wTx-b  with  w=SIW'1(jiIrPi/)  and 
h-=wTb,0  where  b,o=0.5(piyt+pJi).  We  can  also  use  two  output  neurons  to  derive  a 
class-dependent  representation.  Rearrange  the  output  by 


y  =  wTx  -  6  =  wjx  -  6*  -  (wjx  -  &,)=  yt - 


y i 


(6) 
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where  w*=S„/‘(mr  b^),  w^S^'fiw-  bl0),  &*=*ViTbzo,  and  we  have 

>'*=w/r-6*  and_Vr=W(Tx-i/-  Fig.  1  illustrates  such  an  interpretation. 

Once  again,  the  class-dependent  Fisher’s  linear  discriminant  w i  and  related 
bias  bi  derived 
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Figure  3.  Class  separability  visualized  in  the  top  3  PCA  and  wFC- 
DCA  spaces. 


output.  For  example,  for  a  three-class  problem,  it  would  be  straightforward  to  have 
",r=SI1,'l(n,i-bIo),  w2=SIJ(|iI,-b5Q),  w3=SIW'1(pX3-bIo),  6,=w,Tb[0,  62=w2Tbj0,  and 
ij=w3Tb,o.  Fig.  2  illustrates  such  a  case.  Notice  that  such  an  initialization  is  readily 


applicable  to  single  layer  perceptrons. 


EXPERIMENTAL 

VALIDATION 

So  far,  we  have  described  the 
theoretical  observations  on  the 
connections  between  MLP  and 
wFC-DCA,  and  have  presented  the 
hypotheses  of  utilizing  wFC-DCA 
to  optimize  MLP.  We  shall  now 
illustrate  the  effectiveness  of  the 
proposed  methods. 

The  first  experiment  reports 
the  test  on  the  Iris  data  set  (m2= 3 
classes,  50  samples  per  class,  m o=4 
dimensions).  We  compared  the 
class  separability  by  wFC-DCA  to 
that  by  the  popular  principle 
component  analysis  (PCA)  [5].  The 
improvement  was  obtained  by 
15%.  Fig.  3  illustrates  the  class 
separability  in  three-dimensional 


projection  of  gsn*  tubipace  to  top  3  wOCA  span 
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Fisher  10  4211;  wFj»h*rl7.«5l 

Figure  4.  Data  plots  in  the  top  wFC-DCA  and 
MLP  extracted  pseudo-gene  spaces. 


plot.  We  further  used  wFC-DCA  to  initialize  the  MLP  and  compare  the 
classification  results  with  that  of  the  randomly  initialized  MLP.  In  average  across 
wi=2,3,4,  the  classification  accuracy  increases  from  62%  to  66%. 
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We  then  tested  wFC-DCA-MLP  method  on  a  publicly  available  SRBCT 
microarray  gene 
expression  data  set 
consisting  of  w2=  4 
different  small 
round  blue  cell 
tumors  [5].  From 
total  2308 

expressed  genes, 
we  selected  m0=S~8 
genes  for  the 
further  study.  We 
first  used  single 
layer  perceptions 

with  wFC-DCA  extracted  pseudo-genes  as  inputs  and  compared  the  classification 
performance  to  that 
of  using  PCA  as 
previously  reported 
[5].  The  results 
show  that,  using 
only  top  3~5  such 
pseudo-genes, 

100%  prediction 
accuracy  on  the 
total  64  samples 
estimated  by  a  3- 
fold  cross- 

validation  can  be 

achieved.  We  further  tested  a  wFC-DCA  optimized  MLP  and  compared  the 
classification  performance  with  that  of  the  randomly  initialized  MLP  on  the  same 
data  set.  The  results  consistently  show  that,  using  only  3  top  pseudo-genes,  i.e.,  3 
hidden  neurons,  the  classification  accuracy  has  been  increased  from  85%  to  100% 
on  the  total  64  samples  estimated  by  a  3-fotd  cross-validation.  Fig.  4  shows  the 
close  yet  improved  patterns  of  data  plots  in  the  pseudo-genes  spaces  that  were 
initialized  by  wFC-DCA  and  refined  by  the  synaptic  weights  learning  in  the  MLP. 
As  expected,  we  have  also  found  that  these  3  dominant  pseudo-genes  contained 
more  than  80%  of  the  total  intrinsic  Bayes  class  separability  [2]. 

To  obtain  more  quantitatively  comparisons,  we  applied  our  method  to  a 
simulated  data  set  with  mo=4  and  m2=  3.  The  performance  was  estimated  via  a  3- 
fold  cross-validation.  We  compared  the  performance  of  the  MLP  with  random  and 
wFC-DCA  initializations  in  terms  of  the  averaged  misclassification  rate  and  its 
standard  deviation,  with  200  trials  for  each  case.  Table  1  reports  the  validation 
results.  These  controlled  experiments  shown  that  with  wFC-DCA  initialization, 
MLP  reached  faster  and  more  global  convergence  in  that  it  used  much  less  epochs 
and  reached  higher  and  more  robust  prediction  accuracy. 


Misclassification 
rate  (mean) 

Misclassification 
rate  (standard 
deviation) 

Mean  squared 
error  during 
training 

MLP  with 

wFC-DCA 

initialization 

5.23% 

0.0635 

0.0032 

MLP  with 

random 

initialization 

11.2% 

0.1672 

0.0755 

Table  2,  Classification  performance  of  the  MLP  with  wFC- 
DCA  and  random  initializations  on  the  SRBCT  data  set. 


Misclassification 
rate  (mean) 

Misclassification 
rate  (standard 
deviation) 

Mean  squared 
error  during 
training 

MLP  with 

wFC-DCA 

initialization 

0.47% 

0.0141 

0.0169 

MLP  with 

random 

initialization 

1.4% 

0.0656 

0.0279 

Table  1 .  Classification  performance  of  the  MLP  with  wFC- 
DCA  and  random  initializations  on  the  simulated  data  set. 
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As  an  example  of  more  challenging  problems  and  more  rigorous  efforts,  we 
conducted  additional  validations  on  the  SRBCT  data  set  now  with  m»= 20  top 
discriminatory  genes  to  assess  the  impact  of  the  proposed  optimization  on  easing 
the  curse  of  dimensionality.  Once 
again,  the  performance  was 
estimated  via  a  3-fold  cross- 
validation,  and  we  compared  the 
performance  of  the  MLP  with 
random  and  wFC-DCA 
initializations  in  terms  of  the 
averaged  misclassification  rate  and 
its  standard  deviation,  with  200 
trials  for  each  case.  Table  2  reports 
the  validation  results.  From  Fig.  5, 
it  can  be  seen  that  when  the  MLP  is 
initialized  by  the  wFC-DCA  that 
has  an  initial  lower  mean  squared 
error  (MSE)  of  0.00000388,  the 
refined  pseudo-genes  via  the 
converged  weights  of  the  MLP 
were  very  close  to  the  wFC-DCA 
extracted  ones.  Not  surprisingly,  when  the  MLP  is  initialized  randomly  with  a  high 
initial  MSE  of  0.1279,  the  converged  weights  of  the  MLP  were  trapped  into  the 
local  optimum  and  deviated  from  the  correct  ones.  Notice  that  the  correctly 
extracted  discriminatory  pseudo-genes  shall  support  a  higher  separability  between 
the  classes. 

As  a  final  test,  we  considered  the  ALL-AML  dataset  of  the  acute  leukemias 
where  mo=20  and  mi- 4.  Since  the  sample  size  is  relatively  small,  we  used  a  5-fold 
cross  validation  to  estimate  the  classification  performance.  Once  again,  we 
compared  the  performance  of  the  MLP  with  random  and  wFC-DCA  initializations 
in  terms  of  the 
averaged 
misclassification 
rate  and  its  standard 
deviation,  with  200 
trials  for  each  case. 

Table  3  reports  the 
validation  results 
and  Fig.  6  provides 
the  further  evidence 
of  the  positive 
impact  of  the  wFC- 
DCA  initialization 
on  the  performance  of  the  MLP. 

Thus,  we  see  that  wFC-DCA  based  MLP  optimization  works  reasonably  well 
as  expected.  When  the  distributions  of  the  data  set  are  close  to  SFNM  that  can  be 


Misclasstfication 
rate  (mean) 

Misclassification 
rate  (standard 
deviation) 

Mean  squared 
error  during 
training 

MLP  with 

wFC-DCA 

initialization 

3.57% 

0.0301 

0.0684 

MLP  with 

random 

initialization 

22.6% 

0.2446 

0.4862 

Table  3.  Classification  performance  of  the  MLP  with  wFC- 
DCA  and  random  initializations  on  the  ALL-AML  data  set. 


wFC-DCA 
Extracted 
2-D  plot 


MLP 

Extracted 
2-D  plot 


Figure  5.  SRBCT  plots  in  the  top  wFC-DCA 
and  MLP  extracted  pseudo-gene  spaces  (left: 
by  wFC-DCA  initialization.  Right:  by 
random  initialization). 
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effectively  visualized  and  estimated  using  our  Visual  and  Statistical  Data  Analyzer 
(VISDA)  software  (http://www.cbil.ece.vt.edu/software)  {4],  it  gives  excellent 
results.  When  the  distributions  of  the  data  set  are  multimodal  in  nature  with  non- 
Gaussian  classes,  the  optimization  is  less  effective  but  still  quite  reasonable. 
However,  if  the  classes  are  highly  nonlinear  separable,  more  care  shall  be  taken  to 
fully  utilize  the  proposed  method.  While  the  effectiveness  of  this  method  may  be 
data-dependent,  we  would  expect  it  to  be  a  very  useful  tool  in  the  design  of  the 
MLP  classifiers,  given  the  unavailability  of  the  general  case-specific  guidelines 
[6). 


DISCUSSIONS 


We  have  shown  a  close  connection  between  classification  by  MLP  and 
classification  by  LDA.  By  suggesting  a  wFC-DCA  based  optimization  procedure 
for  improving  the  utility  of  MLP,  our  main  contribution  lies  in  the  theoretical 
insight  and  experimental  validation 
on  what  MLP  is  doing. 

One  important  consideration 
with  the  present  MLP  optimization 
is  that  the  complexity  of  the  MLP 
(e.g.,  total  number  of  the  freely 
adjustable  synaptic  weights  in  the 
network)  should  be  minimized  in 
order  to  ease  the  curse  of 
dimensionality  [1],  From  our  later 
experimental  results,  it  can  be  seen 
that  the  estimated  generalizable 
performance,  as  expected, 
decreased  with  the  increase  in  the 
input  dimensions  (e.g.,  m0=20 
versus  m<f=5~8). 

It  must  be  noted,  however,  that 
LDR/LDA  belongs  to  probability 
density  based  method  that  assumes 
Gaussian  model  for  each  class  by  model  estimation,  while  MLP  belongs  to 
geometric  approach  that  finds  the  decision  boundary  directly  from  maximizing 
classification  accuracy.  From  Vapnik’s  philosophy  of  learning  theory  [1,3],  the 
further  trained  MLP  often  outperforms  the  initial  wFC-LDA.  Our  experimental 
results  also  support  such  an  expectation,  particularly  when  the  data  distributions 
deviate  from  SFNM  or  the  sample  size  is  simply  too  small. 

Although  wFC-DCA  is  successful  in  optimizing  the  MLP  in  multiclass 
classification,  this  success  does  not  imply  that  the  method  can  be  effective  for  any 
data  set,  such  as  multiclass  nonlinearly  separable  cases.  A  classification  problem 
could  be  an  intrinsically  nonlinear  problem,  or  becomes  a  nonlinear  problem  after 
dimensionality  reduction  according  to  Cover’s  theorem  on  the  separability  of 


Figure  6.  ALL-AML  plots  in  the  top  wFC- 
DCA  and  MLP  extracted  pseudo-gene  spaces 
(left:  by  wFC-DCA  initialization.  Right:  by 
random  initialization). 
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patterns  [3],  As  aforementioned,  therefore,  the  hidden  layer  of  the  MLP  needs  to 
perform  the  additional  function  of  transforming  a  nonlinearly  separable  problem 
into  a  linear  classification.  This  may  be  achieved  by  the  existing  hidden  layer 
through  dual-purpose  training,  or  one  additional  hidden  layer  shall  be  introduced 
[3],  An  elegant  yet  simple  method  is  to  apply  divide-and-conquer  principle  to  the 
data  set  and  accordingly  introduce  some  pseudo-classes  to  the  output  layer,  such 
that  all  class-pairs  are  linearly  separable  thereafter.  Notice  that  the  discrete 
decision  fusion  can  be  readily  and  effortlessly  done  without  using  any  combiner, 
since  the  pseudo-classes  belong  to  some  of  the  known  classes  as  a  priori.  It  shall 
not  escape  our  notice  that  a  net  reduction  in  MLP  complexity  can  still  be  achieved 
when  mo  is  large,  since  the  total  number  of  weights  in  a  three-layer  MLP  is  mi(mo+ 
m2)  such  that  die  reduction  due  to  m,  surpasses  the  generally  limited  increase  due 
to  m2.  Further  refinements  include  methods  that  allow  a  co-determination  of  m , 
and  m2  for  an  optimum  generalizable  performance  when  the  curse  of 
dimensionality  is  the  issue  [1,3, 4, 5, 7]. 


Appendix 

Computation  of  Stw'm: 

Let  the  eigenvalue  decomposition  of  Sw  be  Stw=EDET  where  E  is  the 
orthogonal  matrix  of  eigenvectors  of  S*  and  D  is  the  diagonal  matrix  of  its 
eigenvalues,  then  Sh/m=ED'l,IET  where  the  matrix  D'1'2  is  computed  by  a  simple 
componentwise  operation  [8], 
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ABSTRACT 

Antiestrogens  induce  both  cytostasis  (cell  cycle  arrest)  and  apoptosis, 
but  the  relationship  between  these  end  points  and  the  signaling  that 
regulates  their  induction  are  unclear.  We  have  previously  implicated 
the  transcription  factor  and  putative  tumor  suppressor  IFN  regulatory 
factor-1  (1RF-1)  in  acquired  antiestrogen  resistance  (Gn  etai,  Cancer  Res, 
62:  3428-3437,  2002).  We  now  show  the  functional  significance  of  1RF-I 
in  affecting  antiestrogen  responsiveness  in  estrogen  receptor-positive 
antiestrogen-sensitive  models  (MCF-7,  T47D,  and  ZR-75-1),  a  model  of 
acquired  antiestrogen  resistance  (MCF7/LCC9;  estrogen  receptor  posi¬ 
tive),  and  a  model  of  de  novo  antiestrogen  resistance  (TV1DA-MB-231; 
estrogen  receptor  negative).  Basal  1RF'-1  mRNA  expression  is  lower  in 
MCF7/LCC9  cells  when  compared  with  IVICF-7,  T47D,  and  ZR-75-1  cells. 
IRF-1  transcriptional  activity  in  MCF-7/LCC9  cells  is  18-fold  lower  than 
that  seen  in  the  parental  cells  (MCF-7/LCC1)  and  is  comparable  with  that 
in  MDA-MB-231  cells.  Although  IRF-1  mRNA  expression  is  induced  by 
ICI  182,780  in  sensitive  cells,  this  regulation  is  lost  in  MCF-7/LCC9  and  is 
absent  in  MDA-MB-231  cells.  Loss  of  IRF-1  regulation  appears  specific  to 
antiestrogen  resistance — resistant  cells  induce  IRF-1  mRNA  in  response 
to  the  cytotoxic  drug  doxorubicin.  A  dominant-negative  IRF-1  eliminates 
the  ICI  182,780-induced  apoptotic  response  (reduced  >4-foid)  and  re¬ 
duces  MCF-7  and  T47D  cell  sensitivity  to  the  antiproliferative  effects  of 
ICI  182,780.  This  effect  is  not  mediated  by  changes  in  cell  cycle  distribu¬ 
tion;  rather,  dominant-negative  IRF-1  reduces  ICI  182,780-induced  apo¬ 
ptosis.  These  data  identify  a  novel  mechanism  of  antiestrogen  resistance 
and  implicate  IRF-1  as  a  key  component  in  signaling  some  ER-mediated 
effects  on  apoptosis/ccll  survival. 

INTRODUCTION 

For  many  women,  antiestrogen  therapy  is  the  least  toxic  and  most 
effective  means  to  manage  their  hormone-dependent  breast  cancer. 
The  most  widely  studied  antiestrogen  has  been  tamoxifen  (TAM), 
which  can  increase  both  disease-free  and  overall  survival  in  breast 
cancer  patients,  reduce  the  incidence  of  estrogen  receptor-positive 
(ER+)  disease  in  high-risk  women,  and  reduce  the  rate  of  bone  loss 
from  osteoporosis  (1,  2).  Although  first  line  antiestrogen  therapy 
remains  the  standard  of  care  for  these  patients  (3-5),  approximately 
one-third  of  all  ER+  breast  tumors  exhibit  de  novo  antiestrogen 
resistance,  and  most  initially  responsive  tumors  eventually  acquire 
resistance  (6). 

The  steroidal  antiestrogen  ICI  182,780  (Faslodex;  Fulvestrant)  has 
successfully  completed  clinical  trials  and  exhibits  considerable  poten¬ 
tial  for  more  widespread  clinical  use  (7).  Most  currently  available 
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antiestrogens  show  little  or  no  significant  activity  in  TAM-resistant 
disease,  which  is  often  treated  with  a  second-line  aromatase  inhibitor. 
However,  ICI  182,780  is  clearly  active  in  patients  who  have  received 
TAM  treatment  and  eventually  recurred  (8).  Furthermore,  two  Phase 
III  clinical  trials  in  TAM-resistant  patients  have  shown  ICI  182,780  to 
be  at  least  as  effective  as  the  potent  aromatase  inhibitor  anastrazole  (9, 
10).  Unlike  most  other  antiestrogens,  ICI  182,780  is  a  pure  ER 
antagonist  (11)  that  can  induce  degradation  of  ER  protein  (12)  and 
inhibit  receptor  dimerization  (13).  Furthermore,  ICI  1 82,780  is  devoid 
of  the  uterotropic  activity  associated  with  the  ability  of  TAM  to 
increase  the  risk  of  developing  endometrial  cancers  (8,  14). 

Antiestrogen  and  estrogen  responsiveness  are  complex  phenotypes, 
and  both  genomic  and  nongenomic  activities  are  functionally  impli¬ 
cated  (6,  15).  In  sensitive  cells,  antiestrogens  are  clearly  cytostatic, 
inducing  a  G0-G,  cell  cycle  arrest  in  vitro.  Clinically,  the  ability  of 
antiestrogens  to  induce  significant  reductions  in  tumor  size  and  in¬ 
creases  in  overall  survival  (1,  2)  and  to  inhibit  the  development  of 
ER+  tumors  in  the  chcmopreventive  setting  (16,  17)  strongly  suggest 
that  these  drags  also  may  be  cytotoxic.  Evidence  implicates  an  induc¬ 
tion  of  apoptotic  cell  death  as  the  major  mechanism  through  which 
antiestrogens  might  induce  a  cytotoxic  effect  (6).  However,  the  rela¬ 
tionship  between  growth  arrest  and  apoptosis  and  how  antiestrogens 
functionally  affect  cell  signaling  to  regulate  these  two  end  points 
remains  to  be  firmly  established. 

It  is  becoming  apparent  that  antiestrogen  resistance  in  ER+  tumors 
is  unlikely  to  be  driven  by  a  single  gene/signaling  pathway.  Thus,  we 
have  invoked  a  gene  network  hypothesis  that  confers  diversity  in 
estrogen/antiestrogen-initiated  signaling  (15,  18,  19).  Ultimately,  we 
envision  multiple  concurrent  signals  through  this  network  of  inte¬ 
grated  and  potentially  interdependent  pathways,  some  antiapoptotic 
and  some  proapoptotic,  with  cellular  response  reflecting  the  dominant 
signals.  In  antiestrogen-unresponsive  cells,  we  hypothesize  that  the 
endocrine  regulation  and/or  function  of  key  components  of  this  net¬ 
work  is  changed  and  that  proapoptotic  signals  are  no  longer  induced 
and  antiapoptotic  signals  have  become  dominant.  To  begin  identifying 
key  genes  that  may  make  up  such  a  network,  we  have  applied  both 
serial  analysis  of  gene  expression  and  gene  expression  microarray 
analyses  to  a  series  of  antiestrogen-sensitive  and  -resistant  cells. 
Among  the  key  genes  identified  is  the  putative  tumor  suppressor  gene 
IFN  regulatory  factor-1  (IRF-1;  Ref.  19). 

Although  initially  identified  as  an  IFN-responsive  gene,  IRF-1  has 
shown  activity  as  a  tumor  suppressor  in  several  studies  (20-22).  For 
example,  IRF-1  is  deleted  in  some  cancers  (23,  24),  and  loss  of  IRF-1 
significantly  increases  tumorigenicity  in  mouse  models  driven  either 
by  ras  or  loss  of  p53  (25).  IRF-1  can  signal  to  apoptosis  in  a 
p53-dependent  or  -independent  manner  (26,  27);  with  or  without 
induction  of  p21  wafl/c,pl  (26)  orp27klpl  (28);  and  through  caspase-1 
(27),  caspase-7  (29),  caspase-8  (30),  and/or  Fas  ligand  (3 1 ).  Loss  of 
p53  activity  is  common  in  breast  cancer  (32).  Nonetheless,  many 
breast  cancers  are  initially  responsive  to  cytotoxic  drags  and  hor¬ 
mones  (1,  33),  implying  that  drug-induced  apoptosis  likely  occurs 
through  both  p53-dependent  and  -independent  mechanisms.  TAM- 
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induced  growth  arrest  can  occur  independently  of  p53  (34),  but  the 
precise  signaling  responsible  for  these  effects  requires  additional 
study.  The  primary  mechanisms  of  cell  growth  arrest  and  apoptosis  for 
ICI  182,780  and  the  importance  of  signaling  through  p53  arc  un¬ 
known. 

In  addition  to  our  previous  study  implicating  IRF-1  in  affecting 
antiestrogen  responsiveness  in  MCF-7  cells  (19),  Harroch  et  at.  (35) 
observed  that  interleukin  6  inhibited  proliferation  and  induced  IRF-1 
mRNA  and  IRF-1  binding  to  its  target  DNA  sequence  in  T47D  cells. 
In  an  immunohislochemistry  study  of  IRF-1  expression  in  breast 
cancer,  the  authors  report  less  IRF-1  expression  in  neoplastic  com¬ 
pared  with  normal  human  breast,  consistent  with  reduced  expression 
of  a  putative  tumor  suppressor  gene.  However,  IRF-1  expression  was 
not  assessed  in  association  with  established  prognostic  markers  or 
clinical  outcome  (36). 

In  this  study,  we  used  the  ER—  MDA-MB-231  cells  as  a  model  of 
de  novo  antiestrogen  resistance  (6).  To  model  de  novo  antiestrogen 
sensitivity,  we  used  the  estrogen-dependent,  ER+  MCF-7  (37)  and 
T47D  cells  (.38)  and  the  ER+,  antiestrogen-sensitive  but  estrogen- 
independent  MCF-7  variant  MCF7/LCC1  (39).  As  a  model  of  ac¬ 
quired  ICI  182,780  resistance,  we  studied  the  MCF7/LCC9  cells, 
which  were  derived  from  MCF7/LCC1  cells  and  are  ER+,  estrogen 
independent,  and  ICI  1 82,780  and  TAM  cross-resistant  (40).  These 
studies  strongly  implicate  signaling  through  IRF-1  and  its  protein 
partners  as  a  critical  mediator  of  antiestrogen  signaling  and  as  a  key 
gene  in  a  broader  gene  network  (15,  19,  41).  Hence,  we  now  show  that 
IRF-1  mRNA  expression  is  induced  by  ICI  182,780  and  repressed  by 
estrogens  in  antiestrogen-sensitive  cells.  Hormonal  regulation  of 
IRF-1  is  absent  in  ER-  cells  and  is  specifically  lost  in  ER+  cells  with 
acquired  antiestrogen  resistance.  Both  MCF-7  and  T47D  breast  cancer 
cells  expressing  a  dominant-negative  IRF-1  (dnlRF-1)  exhibit  a  de¬ 
crease  in  sensitivity  to  ICI  1 82,780.  The  data  separate  the  proapoptotic 
activity  of  ICI  182,780  from  its  ability  to  induce  cell  cycle  arrest  and 
arc  consistent  with  IRF-1  playing  a  critical  role  in  those  proapoptotic 
activities  of  antiestrogens  most  likely  to  contribute  to  their  ability  to 
increase  overall  survival  and  to  reduce  the  risk  of  developing  ER+ 
breast  cancer  (1,  16,  42). 


MATERIALS  AND  METHODS 

Cell  Culture  and  Reagents.  MCF-7,  T47D,  ZR-75-1  (ER+,  estrogen 
dependent,  antiestrogen  sensitive),  and  MDA-MB-231  (ER— ,  estrogen  inde¬ 
pendent,  antiestrogen  unresponsive)  cells  were  routinely  grown  in  improved 
minimal  essential  medium  (IMEM;  Biofluids,  Rockville,  MD)  with  phenol  red 
and  supplemented  with  5%  fetal  bovine  serum.  MCF-7  cells  were  originally 
obtained  from  Dr.  Marvin  Rich  (Michigan  Cancer  Foundation,  Detroit,  MI). 
T47D,  ZR-75-1,  and  MDA-MB-231  cells  were  obtained  from  the  Lombardi 
Comprehensive  Cancer  Center’s  Tissue  Culture  Shared  Resource.  MCF-7/ 
LCC1  (ER+,  estrogen  independent,  antiestrogen  sensitive,  MCF-7  variant; 
Refs.  39  and  43)  and  MCF-7/LCC9  cells  (ER+,  estrogen  independent,  TAM 
and  ICI  182,780  cross-resistant,  MCF-7  variant  derived  directly  from  MCF7/ 
LCC1  by  selection  against  ICI  1 82,780;  Refs.  19  and  40)  were  routinely  grown 
in  IMEM  without  phenol  red  and  supplemented  with  5%  charcoal  stripped  calf 
serum-IMEM.  All  cells  were  maintained  in  a  humidified  incubator  at  37°C  in 
an  atmosphere  containing  95%  air:5%  C02.  The  steroidal  antiestrogen  ICI 
182,780  (Faslodex;  Fulvestrant)  was  kindly  provided  by  Dr.  Alan  Wakeling 
(Zeneca  Pharmaceuticals,  Macclesfield,  Cheshire,  United  Kingdom).  Recom¬ 
binant  human  IFN-y  was  purchased  from  Boehringer  Mannheim  (Mannheim, 
Germany). 

RNA  Extraction.  Total  RNA  was  extracted  from  proliferating  subconflu¬ 
ent  cells  using  the  TRIazol  reagent  (Life  Technologies,  Inc.,  Gaithersburg. 
MD)  according  to  the  manufacturer’s  instructions.  In  brief,  cells  were  rinsed 
with  1  X  PBS  to  remove  serum  and  lysed  by  the  addition  of  the  TRIazol 
reagent.  RNA  was  isolated  by  chloroform  extraction  and  precipitated  using 


isopropanol.  Total  RNA  was  quantified  based  on  the  absorbance  at  260  nM 
using  a  spectrophotometer  (DU640;  Beckman,  Fullerton,  CA). 

Riboprobc  Generation  and  RNase  Protection  Analysis.  The  IRF-1  ribo- 
probe  was  generated  by  reverse  transcriptase-PCR  amplification  of  a  portion  of 
the  IRF-1  mRNA  from  MCF-7  cells.  Amplification  by  PCR  applied  one  95°C 
cycle  for  3  min,  40  cycles  of  95°C  for  30  s,  59°C  for  30  s,  and  72°C  for  1  min 
followed  by  I  cycle  of  72°C  for  5  min,  using  the  following  primers:  forward, 
TCCACCTCTCACCAAGAACC  (bp  533-552),  and  reverse,  TTCCCTTCCT- 
C’ATCCTCATC  (bp  873-892).  To  control  for  equivalent  RNA  loading,  we 
measured  expression  of  the  constitutively  expressed  36B4  mRNA  that  encodes 
the  human  acidic  ribosomal  protein  PO  (44).  The  36B4  riboprobe  was  gener¬ 
ated  as  described  previously  (39). 

RNase  protection  assays  were  performed  as  described  previously  (19,  45). 
In  brief,  plasmids  were  linearized  by  digestion  with  EcoRl  and  transcribed 
with  either  SP6  (IRF-i)  orT7  polymerase  (36B4)  in  the  presence  of  [32P]UTP. 
To  obtain  signals  of  approximately  comparable  intensity,  the  36B4  riboprobes 
were  labeled  with  approximately  one-fifth  of  the  ['12P]UTP  concentration  used 
to  label  the  IRF-1  riboprobes.  The  IRF-1  and  36B4  riboprobes  respectively 
generate  360-  and  220-bp  protected  fragments.  For  each  assay,  30  jug  of  total 
RNA  was  hybridized  to  5  X  I04  dpm  of  probe  for  12-16  h  at  50°C  and 
digested  with  40  p.g/ml  RNase  A  for  30  min  at  25°C.  Digestion  was  terminated 
by  the  addition  of  proteinase  K  (0.25  mg/ml)  and  0.5%  (w/v)  SDS.  Samples 
were  extracted  into  phenobchloroforrmisoamyl  alcohol  (25:24:1)  and  precip¬ 
itated  in  ethanol,  and  the  pellets  were  boiled  in  loading  buffer  and  fractionated 
in  6%  Tris-borate  EDTA-urea  polyacrylamide  gels.  Radioactivity  was  detected 
by  autoradiography  and  quantified  using  Phosphorlmager  analysis  (445SI; 
Molecular  Dynamics,  Sunnyvale,  CA). 

Cell  Lysis  and  Immunoblotting.  For  the  determination  of  IRF-1  protein 
expression,  cells  were  seeded  into  6-well  dishes  at  2  X  105  cells/well  and 
cultured  in  normal  growth  media  for  24  h.  To  examine  the  induction  of  IRF-1 
in  response  to  ICI  182,780  or  estradiol,  cells  were  seeded  at  105  cells/well  1 
day  before  treatment  with  either  1  nM  17j3-estradiol  or  100  nM  ICI  182,780  for 
3  days.  Cells  were  lysed  in  modified  radioimmune  precipitation  assay  buffer 
[150  mM  NaCl,  50  him  Tris,  1%  Igepal  C'A-630,  and  0.5%  deoxycholate  (pH 
7.5)]  supplemented  with  Complete  Mini  protease  inhibitor  mixture  tablets 
(Roche,  Mannheim.  Germany).  Lysates  were  clarified  by  centrifugation,  and 
equal  volumes  were  added  to  2X  Laemmli  sample  buffer  before  boiling  and 
loading  onto  precast  12%  acrylamide  gels  (NuPAGFI  Electrophoresis  System, 
Invitrogen,  Carlsbad,  CA).  Proteins  were  transferred  to  nitrocellulose  mem¬ 
branes  and  incubated  with  primary  antibody  (IRF-1  C-20  at  1 :200;  Santa  Cruz 
Biotechnology,  Santa  Cruz,  CA)  in  Tris  Buffered  Saline  and  Tween-20 
[TBST;  10  mM  Tris  HC1,  150  mM  NaCI,  and  0.05%  Tween-20  (pH  8.0)] 
containing  5%  nonfat  diy  milk  overnight  at  4°C.  Membranes  were  then 
incubated  with  horseradish  peroxidase-conjugated  secondary  antibody  for  1  h 
at  room  temperature  followed  by  enhanced  chemiluminescence  (Amersham 
Biosciences,  Piscataway,  NJ)  and  exposure  to  film  (X-OMAT  Blue  XB-1; 
Kodak,  Rochester,  NY).  To  confirm  equal  loading,  membranes  were  reprolied 
as  described  above  using  a  )3-actin  monoclonal  antibody  (1:5000;  Sigma,  St. 
Louis,  MO). 

Generation  of  dnIRF-1.  Although  small  interfering  RNA  (siRNA)  can  be 
a  powerful  method  to  inhibit  RNAs,  we  did  not  use  this  method  to  block  IRF-1 
because  of  the  recent  reports  of  a  marked  induction  of  an  IFN  response  when 
these  molecules  are  introduced  into  cells  (46,  47).  Thus,  the  best  remaining 
approach  was  the  use  of  a  stably  expressed  dominant-negative  strategy.  A 
wild-type  IRF-1  cDNA  (kindly  provided  by  Dr.  Taniguchi,  University  of 
Tokyo,  Japan)  was  subcloned  into  the  Xho 1  site  of  the  pGEM7Z  expression 
vector  (Promega.  Madison,  WI)  and  linearized  with  Bgl II.  The  dnIRF-1  com¬ 
prises  the  full-length  wild-type  IRF-1  cDNA  with  a  deletion  of  bp  647-1173 
and  contains  both  the  3'  and  5'  untranslated  regions,  the  DNA-binding  domain, 
and  the  nuclear  localization  sequences  of  IRF-1.  We  constructed  dnlRF-1  by 
PCR  amplification  using  bp  630-647  (TAGCAGGGCCCCTGG)  as  the  for¬ 
ward  primers  and  bp  1173-1187  (ATCAGAGAAGGTATCAGG)  as  the  re¬ 
verse  primers.  PCR  conditions  were  30  cycles  of  95°C  for  30  s,  45°C  for  45  s, 
and  72°C  for  5  min.  Integrity  of  the  dnIRF-1  sequence  was  confirmed  by 
standard  dideoxy-mediated  chain-termination  sequencing.  dnIRF-1  was  sub¬ 
cloned  into  both  the  Xho\  site  of  the  pcDNA3  mammalian  expression  vector 
(Invitrogen,  Carlsbad,  CA)  and  into  the  Xho\  site  of  the  pbiEGFP-tet  vector  [a 
plasmid  expressing  enhanced  green  fluorescent  protein  (EGFP)]  under  the 
control  of  a  bidirectional  tetracycline-responsive  promoter  (Clontech,  Palo 
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Alto,  CA),  The  IRF-1  riboprobe  (above)  identifies  dnIRF-1,  generating  a 
1 15-bp  protected  dnIRF-1  fragment. 

Transient  Transfections  and  Luciferase  Reporter  Assay.  Cells  were 

transfected  using  the  FuGENE  6  method  (Roche  Diagnostics,  Indianapolis, 
IN).  For  reporter  assays,  8  x  I04  cells/well  were  plated  in  12-well  plates  and 
allowed  to  grow  for  24  h  before  transfection.  Cells  were  cotransfected  with  the 
pISRE-luc  plasmid  (contains  five  copies  of  the  IFN-stimulated  response  ele¬ 
ment;  ISRE)  as  provided  in  the  PathDetect  kit  (Promega)  and  with  either  a 
standard  control  comprising  a  pcDNA3  plasmid  without  the  ISRE  or  a 
pcDNA3  plasmid  containing  the  cDNA  from  either  IRF-1  or  dnIRF-1.  To 
control  for  transfection  efficiency,  a  pRL-SV40  plasmid  (Promega)  containing 
the  Renilla  luciferase  gene  under  the  control  of  a  constitutive  SV40  promoter 
also  was  cotransfected  into  cells.  One  fxg  of  plasmid  DNA  was  added  to 
serum-free  media  containing  the  FuGENE  6  reagent  and  allowed  to  incubate 
for  30  min  at  room  temperature.  Where  appropriate,  cells  were  maintained  in 
growth  media  with  or  without  500  lU/ml  IFN-y  for  24-48  h.  Subsequently, 
cells  were  lysed,  and  activation  of  the  ISRE-luciferase  construct  was  measured 
using  the  Dual  Luciferase  Assay  kit  (Promega)  according  to  the  manufacturer’s 
instructions.  Luminescence  was  quantified  using  a  Luniat  LB  9501  luminom- 
eter  (EG&G  Berthold,  Bundoora,  Victoria,  Australia). 

Stable  Transfection  with  dnIRF-i.  For  stable  transfections,  cells  were 
plated  in  T-75  cm3  plastic  tissue  culture  flasks  at  a  density  of  0.5  X  I06 
cells/flask  and  grown  for  24  h  before  transfection.  A  total  of  8  p.g  of  plasmid 
DNA  were  transfected  into  MCF-7  cells  stably  transfected  with  the  tetR 
protein  (MCF-7/VPI6;  generously  provided  by  Dr.  Susan  Conrad,  Michigan 
State  University,  East  Lansing,  MI)  orT47D  cells  using  the  FuGENE6  method 
(above).  Cells  were  transfected  with  either  an  empty  pBI-EGFP-tet  plasmid 
containing  the  EGFP  selectable  marker  (Clontech)  or  one  containing  the 
dnIRF-1  cDNA  and  the  pBABF.  plasmid  (kindly  provided  by  Dr.  Matthew 
Ellis,  Washington  University,  St.  Louis,  MO)  encoding  for  puromycin  resist¬ 
ance.  Stably  transfected  cells  were  selected  for  growth  in  the  presence  of  I 
ix g/ml  puromycin.  Puromycin-resistant  colonies  expressing  EGFP  as  measured 
by  standard  fluorescence  activated  cell  sorting  (FACS)  were  expanded  and 
screened  for  expression  of  the  dnIRF-1  by  RNase  protection.  All  of  the 
transfectanls  used  in  these  experiments  were  from  pooled  populations.  Cells 
transfected  with  the  empty  control  vector  were  designated  MCF7/ctrl  and 
T47D/ctrl  and  those  transfected  with  the  dnIRF-1  were  designated  MCF7/ 
dnlRF-1  and  T47D/dnIRF-l. 

Cell  Proliferation.  MCF-7/ctrl.  T47D/ctrl,  MCF7/dnlRF-l,  and  T47D/ 
dnIRF-1  cells  were  sorted  aseptically  by  FACS  in  the  Lombardi  Comprehen¬ 
sive  Cancer  Center  Flow  Cytometry  Shared  Resource  for  EGFP  expression  and 
plated  into  1 2-well  plastic  tissue  culture  plates  at  a  density  of  1 .5  X  I03  MCF-7 
or  2  X  103  T47D  cells/well.  Twenty-four  h  post  plating,  cells  were  treated  with 
100  nM  or  1  juM  ICI  182,780  or  vehicle  control  for  72  h.  Cells  were  then 
trypsinized,  resuspended  in  PBS,  and  counted  in  a  Beckman  Coulter  counter 
(Beckman  Coulter  Corp.,  Fullerton.  CA)  to  assess  cellular  proliferation. 

Cell  Cycle  Analyses.  Cells  stably  transfected  with  the  dnIRF-1  or  empty 
control  plasmids  were  plated  in  T-75  cm2  plastic  tissue  culture  flasks  at  a 
concentration  of  0.5  X  106  and  allowed  to  grow  for  3  days.  Cells  were  then 
analyzed  for  alterations  in  cell  cycle  via  FACS.  FACS  analysis  was  conducted 
by  the  Lombardi  Comprehensive  Cancer  Center  Flow  Cytometry  Shared 
Resource,  according  to  the  method  of  Vindelov  el  al.  (48). 

Apoptosis.  Staining  for  annexin  V,  an  optimal  assay  for  detecting  apoptosis 
in  MCF-7  cells  (49,  50),  was  performed  according  to  the  manufacturer’s 
instructions  (R&D  Systems,  Minneapolis,  MN).  Control-  or  dnIRF-1-trans- 
fected  MCF-7  cells  (1  X  106)  were  seeded  in  T-75  cm2  plastic  tissue  culture 
dishes  and  allowed  to  grow  for  24  h.  Cells  were  then  treated  with  ICI  182,780 
or  vehicle  for  72  h,  trypsinized,  and  pelleted  by  centrifugation.  Cell  pellets 
were  rinsed  twice  in  ice-cold  PBS  and  stained  with  5  p.g/ml  7-aminoactino- 
mycin  D  for  15  min  at  room  temperature.  After  staining  with  7-aminoactino- 
mycin  D,  cells  were  washed  and  resuspended  in  annexin  V  binding  buffer  [10 
mu  4-(2-hydroxyethyl)-l-piperazineethanesulfonic  acid/NaOH,  140  ntM  NaCI, 
and  2.5  mM  CaCI2  (pH  7.4)].  Cells  (1  X  105)  were  then  stained  with  0.3  pg  of 
phycoerythrin-conjugated  annexin  V  in  the  dark,  and  flow  cytometric  analysis 
was  performed  using  a  FACStarplus  flow  cytometer  (Becton-Dickinson,  Moun¬ 
tain  View,  CA)  to  determine  the  proportion  of  apoptotic  cells  in  each  sample. 
Apoptosis  was  measured  only  in  cells  expressing  the  dnlRF-1  transgene  or 
empty  vector  (as  assessed  by  concurrent  EGFP  expression). 


Statistical  Methods.  Student’s  /  test  was  used  to  compare  two  groups  in 
which  the  data  are  normally  distributed;  a  Wilcoxon  /  test  was  used  to  compare 
groups  in  which  data  are  not  normally  distributed.  For  multiple  comparisons, 
ANOVA  was  used  with  a  post  hoc  I  test  for  multiple  comparisons.  Where 
several  experimental  groups  were  compared  with  the  same  control,  we  used 
Dunnett's  test  (51). 


RESULTS 

IRF-1  Is  Differentially  Expressed  in  Antiestrogen-Sensitive  and 
-Resistant  Breast  Cancer  Cells.  We  measured  basal  expression  of 
IRF-1  ntRNA  by  RNase  protection  analyses  in  two  models  of  anties¬ 
trogen  resistance:  MCF-7/LCC9  cells  (acquired  resistance,  ER+,  ICI 
182,780  and  TAM  cross-resistant)  and  MDA-MB-231  (de  novo  re¬ 
sistance,  ER— ,  ICI  182,780  and  TAM  cross-resistant;  Ref.  6).  MCF- 
7/LCC9  cells  exhibit  a  4.2-fold  (P  <  0.05)  and  2. 6-fold  (P  <  0.01) 
lower  expression  of  IRF-1  mRNA  than  the  antiestrogen-sensitive 
MCF-7  and  MCF-7/LCC1  cells,  respectively  (Fig.  I,  A  and  B).  IRF-1 
mRNA  expression  in  MCF-7/LCC9  cells  is  not  significantly  different 
from  that  in  ER—  MDA-MB-231  cells. 

The  half-life  of  the  IRF-1  protein  is  less  than  30  min,  and  changes 
in  mRNA  levels  appear  closely  associated  both  with  changes  in  IRF-1 
protein  expression  and  IRF-1  transcriptional  activity  as  measured 
using  an  ISRE-based  promoter-reporter  assay  (52).  To  confirm  this 
association,  we  performed  immunoblotting  for  the  cells  lines  shown  in 
Fig.  1  ,A  and  B.  The  data  in  Fig.  1C  show  that,  as  expected,  the  mRNA 
and  protein  levels  are  comparable.  The  MCF-7/LCC1  and  MCF-7/ 
LCC9  cells  compared  here  are  derived  from  the  same  parental  MCF-7 
cell  line  and  have  similar  transfection  efficiencies  as  assessed  by 
transfection  with  /3-galactosidase  (not  shown).  Activity  of  the  cotrans¬ 
fected  Renilla  construct  (constitutively  active)  was  used  to  correct  for 
any  minor  differences  in  transfection  efficiency.  MCF-7/LCC9  cells 
exhibit  18-fold  lower  basal  ISRE  activity  than  their  immediate  paren¬ 
tal  cells  MCF-7/LCC1  (P  <  0.01;  Fig.  ID).  These  data  reflect  the 
differential  IRF-1  mRNA  expression  that  we  have  previously  detected 
in  gene  expression  microarrays  (19)  and  now'  confirm  by  RNase 
protection  analyses  (Fig.  1  ,A  and  B)  and  immunoblot  (Fig.  1C).  These 
data  also  confirm  that  IRF-1  mRNA  expression,  protein  expression, 
and  transcriptional  activation  are  closely  related  in  breast  cancer  cells. 
Furthermore,  these  data  show  that  IRF-1  expression  and  activity  are 
significantly  lower  in  ER+  and  ER—  antiestrogen-resistant  cells 
compared  with  antiestrogen-sensitive  cells. 

IRF-1  mRNA  Expression  Is  Regulated  through  ER  in  Breast 
Cancer  Cells.  Data  on  basal  expression  in  the  sensitive  and  resistant 
cells  imply  that  estrogens  and  antiestrogens  may  affect  IRF-1  mRNA 
expression.  Fig.  2,  A  and  B,  shows  the  ability  of  100  nM  ICI  182,780, 
a  clinically  relevant  concentration  (53),  to  induce  IRF-1  mRNA  in  the 
three  best  characterized  and  most  widely  used  ER+  human  breast 
cancer  cell  lines  (6);  MCF-7  (P  <  0.001);  T47D  (P  <  0.001);  and 
ZR-75-1  (P  <  0.05).  This  effect  appears  to  be  ER  mediated,  because 
IRF-1  mRNA  induction  by  ICI  182,780  is  blocked  in  the  presence  of 
estradiol  in  MCF-7  cells  (Fig.  2,  C  and  D)  and  ICI  1 82,780  is  unable 
to  induce  IRF-1  in  the  ER—  MDA-MB-231  cells  (Fig.  2,  A  and  B). 

The  dose-response  relationships  for  the  endocrine  regulation  of 
IRF-1  mRNA  in  MCF-7  and  T47D  cells  were  determined.  Cells  were 
plated  and  24  h  later  treated  w'ith  various  doses  of  ICI  1 82,780  or 
0.1%  (v/v)  ethanol  vehicle  for  72-96  h  before  RNA  isolation.  ICI 
182,780  induces  IRF-1  mRNA  in  a  dose-dependent  manner,  with  a 
maximal  2.5-fold  induction  at  a  dose  of  100  nM  in  both  MCF-7  (Fig. 
3,  A  and  B;  P  <  0.001)  and  T47D  (Fig.  3,  C  and  D;  P  <  0.01)  cells. 
In  contrast,  IRF-1  mRNA  expression  is  significantly  down-regulated 
in  response  to  treatment  with  1  n.M  estradiol  (Fig.  3,  E  and  F; 


P  s  0.001). 
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Fig.  1 .  Basal  IRF-I  mRNA  and  protein  expression  and  transcriptional  activation  in  breast  cancer  cell  lines.  A.  representative  RNase  protection  assay.  36B4.  loading  control.  8,  IRF-I 
mRNA  expression  measured  by  RNase  protection  and  presented  as  mean  ±  SE  of  three  determinations,  in  which  intensity'  is  expressed  as  a  ratio  of  IRF-I  :36B4  (ANOVA,  P  -  0.009; 
*,  P  <  0.05  MCF-7  versus  MCF7/LCC9;  P  <  0.01  MCP-7/I.CCI  versus  MCF7/LCC9;  P  =  0.52S  MDA-MB-23 1  versus  MCF7/I.CC9).  C,  representative  immunoblot of  IRF- 1  protein. 
j3 -actin,  loading  control,  D,  basal  transcriptional  activity  of  IRF- 1  in  breast  cancer  cell  lines  (ISRE-luc  promoter-reporter  assay).  Data  represent  mean  ±  SE  of  four  determinations  and 
are  presented  as  relative  light  units.  ANOVA.  P  <  0.01 ;  *,  P  <  0.001  for  MCF-7  versus  each  of  MCT7ICC I ,  MCF7/LCC9,  and  MDA-MB-23 1 .  *,  P  —  NS  for  MCF7/LCC9  versus 
MDA-MB-23 1. 


ER-Mediated  Regulation  of  IRF-1  Is  Specifically  Lost  in 
Anti-estrogen-Resistant  Cells.  Although  the  data  in  Fig.  2  and  Fig. 
3  show  the  regulation  of  IRF-1  mRNA  expression,  this  would  likely 
be  of  limited  functional  relevance  if  similar  patterns  of  regulation 
occur  in  antiestrogen-resistant  cells.  However,  the  ability  of  both 
estrogen  and  ICI  182,780  to  regulate  IRF-1  mRNA  expression  is  lost 
in  the  MCF-7/LCC9  cells  (Fig.  4,  A-D),  and  ICI  182,780  is  unable  to 
regulate  IRF-1  expression  in  the  antiestrogen-resistant  ER-  ceil  line, 
MDA-MB-23 1  (Fig.  2,  A  and  B).  In  contrast,  the  estrogen-indepen- 
dent  but  antiestrogen-sensitive  MCF7/LCC1  cells  retain  the  ability  of 
estrogen  to  inhibit  IRF-1  tnRNA  expression  (Fig.  4,  A  and  B).  Con¬ 
sistent  with  the  data  in  Fig.  1,  we  found  similar  changes  in  the 
hormonal  regulation  of  IRF-1  protein  in  immunoblots  (not  shown). 
MCF7/LCC1  cells  do  not  require  estrogens  to  grow  in  vitro  or  in  vivo 
(estrogen  independent)  but  retain  some  estrogen  responsiveness  (40, 
43).  Thus,  the  apparent  loss  of  ER-mediated  regulation  of  IRF-1  is 
associated  with  acquired  antiestrogen  resistance  but  not  estrogen 
independence. 

We  then  asked  whether  the  specificity  of  endocrine  regulation  of 
IRF-1  that  is  lost  in  antiestrogen-resistant  cells  also  extended  to 
nonhormonal  inducers  of  IRF-1.  In  mouse  embryonic  fibroblasts, 
IRF-1  is  induced  by  treatment  with  the  cytotoxic  drug  doxorubicin 
(26),  one  of  the  most  active  single  cytotoxic  drugs  used  in  the 
treatment  of  breast  cancer  (54,  55).  MCF-7  cells  induce  IRF-1  mRNA 
in  response  to  treatment  with  doxorubicin  (Fig.  5,  A  and  B).  Impor¬ 
tantly,  MCF-7/LCC9  cells  retain  their  ability  to  regulate  IRF- 1  ex¬ 
pression  in  response  to  1  p,M  doxorubicin  (increase  by  over  7-fold; 
P  £  0.001),  a  response  also  shared  by  MDA-MB-231  cells  (increase 
by  over  4-fold;  P  <  0.05;  Fig.  6,  A  and  B).  Thus,  antiestrogen 
resistance  is  associated  with  a  specific,  ER-mediated  change  in  the 
regulation  of  IRF-1  expression,  rather  than  a  global  loss  of  IRF-1 
mRNA  regulation.  This  specificity  allows  cells  to  retain  the  ability  to 
induce  IRF-1  and  undergo  an  IRF-1 -regulated  apoptotic  cell  death  in 
response  to  other  cytotoxic  agents. 


ICI  182,780-Induced  Inhibition  of  Cell  Proliferation  Is  Reduced 
by  dnIRF-1.  To  study  the  functional  relevance  of  changes  in  IRF-1 
activity  in  affecting  antiestrogen  responsiveness,  we  created  dnlRF-1. 
dnlRF-1  contains  the  IRF-1  DNA  binding  domain  and  nuclear  local¬ 
ization  signal  but  lacks  the  protein  binding  and  transcriptional  acti¬ 
vation  domains.  Dominant-negative  activity  of  dnlRF-1  was  apparent 
in  its  ability  to  inhibit  basal  ISRE  activity,  activity  induced  by  tran¬ 
sient  expression  of  wild-type  IRF-1,  and  IFN-y  stimulated  ISRE 
activity  in  MCF-7  (Fig.  7;  all  comparisons  P  <  0.001;  Student’s  t 
test).  Control  cells  transfected  with  an  empty  expression  vector  exhibit 
no  regulation  of  ISRE  activity;  control  cells  transfected  with  empty 
luc  vector  (no  ISRE)  showed  no  activity  (not  shown).  Importantly,  we 
did  not  want  dnIRF-1  to  eliminate  all  IRF-1  transcriptional  activity 
but  rather  to  inhibit  activity  to  an  extent  broadly  equivalent  to  that 
induced  by  antiestrogens.  Complete  loss  of  IRF-1  could  induce  com¬ 
pensatory  responses  of  uncertain  biological  relevance,  and  there  is  no 
evidence  that  IRF-1  expression  is  fully  lost  in  breast  tumors  (36). 
Because  IRF-1  mRNA  expression  levels  are  closely  related  to  ISRE 
activity  in  breast  cells  (Fig.  1),  the  data  in  Fig.  7  imply  that  dnIRF-1 
inhibits  a  level  of  ISRE  activity  broadly  comparable  with  that  induced 
by  antiestrogens  in  sensitive  cells  (Fig.  2;  Fig.  3). 

dnIRF-1  was  stably  transfected  into  both  MCF-7  and  T47D  cells, 
and  the  ability  of  dnlRF-1  to  affect  ICI  1 82,780-induced  inhibition  of 
cell  proliferation  was  measured  (Fig.  8,  A  and  B).  In  cell  proliferation 
assays,  expression  of  dnIRF-1  significantly  reduced  responsiveness  of 
the  MCF-7/dnIRF-l  and  T47D/dnIRF-l  transfectants  to  ICI  182,780. 
At  100  nM,  the  dose  that  maximally  induces  IRF-1  mRNA  expression, 
MCF-7/dnlRF-l  and  T47D/dnIRF-l  transfectants  were  significantly 
less  sensitive  to  ICI  1 82,780  compared  with  their  respective  controls. 
Similar  differences  in  the  responsiveness  of  MCF-7/dnIRF-l  and 
T47D/dnIRF-l  cells  were  seen  at  1  pM  ICI  182,780  (P  =  0.002, 
MCF-7  100  nM  ICI  182,780;  P  =  0.013,  MCF-7  1  pM  ICI  182,780; 
P  =  0.002,  T47D  100  nM  ICI  182,780;  P  =  0.043,  T47D  1  ptM  ICI 
182,780;  Student’s  t  test;  Fig.  8,  A  and  B). 
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Fig.  2.  Endocrine  regulation  of  IRF-1  mRNA  transcription  in  breast  cancer  coll  lines.  A.  representative  RNase  protection  data  for  IRF-1  mRNA  expression.  B.  ICI  182,780  induces 
IRF-1  mRNA  expression  in  ERr  but  not  ER—  breast  cancer  cells.  Cells  were  treated  with  100  nst  ICI  182,780  or  ethanol  vehicle.  Data  represent  mean  ±  SE  of  three  determinations, 
in  which  absorbance  is  expressed  as  a  ratio  of  IRF-1 :36B4  and  represented  as  a  percentage  of  vehicle  control-treated  cells.  *,  P  —  0.001  for  ICI  182,780  versus  control-treated  cells 
within  each  cell  line  for  MCF-7  andT47Dand  P  <  0.05  forZR75.  P  ~  0.880  for  MDA-MB-231;  Student’s  nest.  C,  representative  RNase  protection  data  for  IRF-1  mRNA  expression; 
D.  estradiol  reverses  ICI  182,780-induced  IRF-1  mRNA  transcription  in  MCF-7  cells.  *,  P  <  0.001;  Student’s  /  test. 


dnIRF-1  Does  Not  Affect  ICI  182,780-induced  Changes  in  Cell 
Cycle  Distribution.  Although  antiestrogens  can  affect  both  ceil  cycle 
distribution  and  the  rate  of  apoptosis,  cell  proliferation  assays  measure 
the  sum  of  these  activities.  Hence,  we  asked  directly  whether  the 
effects  of  dnlRF-1  on  proliferation  reflected  an  inhibition  of  the 
ability  of  1C1  182,780  to  arrest  cells  in  G0-G,.  The  data  in  Fig.  9  show 
that  dnIRF-1  does  not  affect  the  ICI  1 82,780-induced  cell  cycle  arrest 
in  G„-G,  in  either  MCF-7/dnIRF-l  (Fig.  9 A)  or  T47D/dnIRF-l  cells 
(Fig.  9 B).  Thus,  the  residual  antiproliferative  effects  of  ICI  1 82,780 
(Fig.  8)  in  dnlRF-1 -expressing  cells  are  those  conferred  by  cell  cycle 
arrest.  These  data  strongly  implicate  changes  in  apoptosis  as  being  the 
primary  mechanism  through  which  dnIRF-1  reduces  the  antiprolifera¬ 
tive  effects  of  ICI  182,780  in  MCF-7  and  T47D  cells. 

ICI  182,780-induced  Apoptosis  Is  Reduced  by  dnIRF-1.  To 
determine  whether  the  effects  of  dnIRF-1  on  cellular  sensitivity  to  ICI 
182,780  are  mediated  by  its  ability  to  influence  signaling  to  apoptosis, 
the  ability  of  dnIRF-1  to  affect  ICI  182,780-induced  apoptosis  was 
assessed  directly  by  measuring  annexin  V  staining  (49,  50).  Apoptosis 
was  measured  only  in  those  cells  expressing  the  dnIRF-1  transgene  or 
empty  vector  control  (as  assessed  by  EGFP  expression),  to  ensure  that 
any  effects  were  likely  to  be  a  direct  result  of  the  inhibition  of  IRF-1. 
When  treated  with  100  nM  ICI  182,780,  30%  of  MCF-7  control 
transfectants  undergo  apoptosis.  Expression  of  dnIRF-1  significantly 
reduces  this  ICI  182,780-induced  apoptotic  response  by  >4-fold  in 
MCF-7/dnIRF-l  cells  to  7%  (Fig.  10;  P  <  0.034).  The  basal  rate  of 
apoptosis,  measured  in  control  MCF-7  cells  treated  with  ethanol,  is 
about  5%  (Fig.  10).  Thus,  the  full  apoptotic  response  to  ICI  182,780 
is  blocked  by  dnIRF-1.  Studies  were  also  performed  in  T47D  and 


T47D/dnIRF-l  cells,  which  unlike  MCF-7  cells  contain  a  mutant  and 
nonfunctional  p53.  T47D/dnIRF-l  cells  show  a  similar  4-fold  reduc¬ 
tion  in  the  ability  of  ICI  182,780  to  induce  apoptosis  (not  shown). 
These  data  with  dnIRF-1  in  both  MCF-7  (wild-type  p53)  and  T47D 
(mutant  p53)  strongly  suggesting  that  IRF-1  is  a  critical  mediator  of 
this  signal  and  that  its  activities  are  independent  of  p53. 

DISCUSSION 

Data  from  our  previous  studies  demonstrate  an  association  between 
IRF-1  expression  and  acquired  cross-resistance  to  antiestrogens  (19). 
We  now  show  that  IRF-1  is  a  key  signaling  protein  involved  in 
mediating  the  sensitivity  of  breast  cancer  cells  to  ICI  1 82,780-induced 
apoptosis.  Basal  IRF-1  mRNA  expression  is  down-regulated  in 
antiestrogen-resistant  cells  (both  the  ER+  MCF7/LCC9  model  of 
acquired  resistance  and  the  ER—  MDA-MB-231  model  of  de  novo 
resistance).  Moreover,  the  ability  of  antiestrogens  to  regulate  IRF-1 
mRNA  transcription  is  absent  in  the  ER-  and  lost  in  the  ER+ 
antiestrogen-resistant  cells.  The  functional  relevance  of  these  obser¬ 
vations  is  shown  by  the  ability  of  dnlRF-1  to  reduce  significantly  the 
antiproliferative  effects  of  ICI  182,780  in  both  the  antiestrogen- 
sensitive  MCF-7  and  T47D  cells.  Notably,  dnIRF-1  does  not  eliminate 
basal  IRF-1  activity  in  these  cells.  Loss  of  IRF-1  activity  could  induce 
confounding  compensatory  effects  unlikely  to  occur  in  breast  tumors, 
which  appear  to  retain  detectable  IRF-1  protein  expression  (36).  Thus, 
the  data  reported  herein  likely  reflect  the  contribution  of  only  the 
antiestrogen  induced  IRF-1. 

Specificity  of  these  effects,  in  the  context  of  the  signaling  of  IRF-1 
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Fig.  3.  Dose-response  and  time-response  relationships  for  the  induction  of  IRF-I  mRNA  by  ICI  182,780  in  MCF-7  ( A  and  B)  and  T47D  cells  (C  and  D).  A,  representative  RNase 
protection  assay.  36B4,  loading  control.  B,  dose-response  relationship  for  IRF-1  mRNA  induction  in  MCF-7  cells.  Cells  were  treated  with  ICI  182,780  or  ethanol  vehicle  for  72  h.  Data 
represent  mean  ±  SE  of  three  independent  replicate  experiments,  in  which  absorbance  is  expressed  as  a  ratio  of  IRF-1 :36B4.  *,  /»  <  0.01  for  treatments  versus  control;  Dunnett’s  test. 
C,  representative  RNase  protection  assay.  36B4,  loading  control.  D,  dose-response  relationship  for  IRF-1  mRNA  induction  in  T47D  cells.  Cells  were  treated  as  in  B.  Data  represent 
mean  ±  SE  of  three  independent  replicate  experiments,  in  which  absorbance  is  expressed  as  a  ratio  of  IRF-1 :36B4.  *,  P  ^  0.01  for  treatments  versus  control;  Dunnett’s  test.  E, 
representative  RNase  protection  assay.  36B4.  loading  control.  F,  estradiol  regulation  ofIRF-1  mRNA  expression  in  MCF-7  cells.  Cells  were  stripped  of  estrogens,  grown  in  the  absence 
of  estrogen  (charcoal  stripped  calf  senim-IMEM),  and  then  treated  with  either  1  nM  estradiol  or  ethanol  vehicle  for  24  h.  Data  represent  mean  ±  SE  of  three  independent  replicate 
experiments,  in  which  absorbance  is  expressed  as  a  ratio  of  IRF-1 :36B4  and  represented  as  a  percentage  of  vehicle  control-treated  cells.  *,  P  <  0.001;  Student’s  t  test. 


in  response  to  ER-mediated  events,  also  is  apparent.  Estradiol  blocks 
antiestrogen-induced  IRF-1  in  MCF-7  cells,  and  no  endocrine  regu¬ 
lation  is  seen  in  ER—  cells.  Neither  the  loss  of  endocrine  regulation  in 
MCF7/LCC9  cells  nor  the  absence  of  its  endocrine  regulation  in 
MDA-MB-231  cells  compromises  the  ability  of  doxorubicin  to  induce 


IRF-1  in  these  cells  and  to  inhibit  their  proliferation  (not  shown).  This 
latter  effect  is  consistent  with  patterns  of  clinical  responses  to  these 
drugs,  because  breast  cancer  patients  who  are  resistant  to  antiestro¬ 
gens  can  respond  to  cytotoxic  chemotherapy. 

The  ability  of  dnIRF-1  to  block  ICI  182,780-induced  inhibition  of 
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Fig.  4.  Hormonal  regulation  of  IRF-1  expression  is  lost  in  antiestrogen-resistant  cells.  A,  representative  RNasc  protection  assay.  36B4,  loading  control.  B,  estradiol  regulation  of  IRF-I 
mRNA  expression.  Cells  were  stripped  of  estrogens,  grown  in  the  absence  of  estrogen  (charcoal  stripped  calf  serum-IMEM),  and  then  treated  with  either  1  nM  estradiol  or  ethanol  vehicle 
for  24  h.  Data  represent  mean  ±  SE  of  three  independent  replicate  experiments,  in  which  absorbance  is  expressed  as  a  ratio  of  IRF-I  ;36B4  and  represented  as  a  percentage  of  vehicle 
control -treated  cells.  *,  P  <  0.001  for  estradiol  versus  control  for  each  of  MCF7  and  MCF7/LCC1  cells:  P  —  0.864  for  estradiol  versus  control-treated  MCT-7/LCC9  cells;  Student’s 
t  test.  C,  ICI  182,780  regulation  of  IRF-1  mRNA  expression;  representative  RNase  protection  assay.  36B4,  loading  control.  Cells  were  treated  with  100  n\i  ICI  182,780  or  ethanol  vehicle 
for  72  h.  D,  ICI  182,780  regulation  of  IRF-1  mRNA  expression.  Data  represent  mean  ±  SE  of  three  determinations,  in  which  absorbance  is  expressed  as  a  ratio  of  IRF-1 :36B4  and 
represented  as  a  percentage  of  vehicle  control-treated  cells.  *,  P  ^  0.001  for  ICI  182,780  versus  control-treated  MCF-7  cells;  P  =  0.971  for  ICI  182,780  versus  control-treated 
MCF7/LCC9  cells;  Student’s  /  test. 


cell  proliferation  in  MCF-7/dnlRF-l  and  T47D/dnIRF-l  cells  could 
reflect  changes  in  the  effects  of  ICI  182,780  on  cell  cycle  and/or 
apoptosis.  However,  dnlRF- 1  does  not  affect  ICi  182,780-induced  cell 
cycle  arrest  when  expressed  in  either  MCF-7  or  T47D  cells.  In  marked 
contrast,  dnIRF-1  effectively  eliminates  ICI  182,780-induced  apopto¬ 
sis  in  both  MCF-7/dnlRF-l  and  T47D/dnIRF-l  cells.  Thus,  we  can 
separate  cell  cycle  arrest  from  apoptosis  and  attribute  a  significant 
component  of  antiestrogen-induccd  apoptotic  signaling  to  IRF-1.  Be¬ 
cause  dnIRF-1  abrogates  ICI  182,780-induced  apoptosis  (Fig.  10) 
while  enhancing  cell  growth  by  approximately  50%  (Fig.  8),  apoptosis 
and  cell  cycle  arrest  likely  contribute  equally  to  the  apparent  antipro¬ 
liferative  effects  of  ICI  182,780. 

Functionally  separating  antiestrogen-induced  apoptosis  from  anti¬ 
estrogen-induced  growth  arrest  has  several  important  implications. 
Novel  therapeutic  approaches  designed  to  increase  the  proapoptotic 
effects  of  antiestrogens  may  be  an  effective  means  to  improve  their 
ability  to  increase  overall  survival  in  patients  because  this  should 
increase  the  proportion  of  cells  undergoing  apoptotic  cell  death. 
Modalities  that  increase  only  the  ability  of  anliestrogens  to  induce  a 
cell  cycle  arrest  will  likely  be  a  less  effective  strategy.  For  example, 
many  arrested  cells  will  survive  and  thereby  have  more  opportunities 
to  adapt,  acquire  resistance,  and  generate  subsequent  disease  recur¬ 
rence.  Measuring  basal  IRF-1  expression  and/or  the  ability  of  an 
antiestrogen  to  induce  IRF-1  in  the  neoadjuvant  setting  may  improve 
the  prediction  of  endocrine  responsiveness.  Currently,  we  incorrectly 
predict  antiestrogen  sensitivity  in  66%  of  ER+/progesterone  receptor¬ 
negative,  55%  of  ER-/progesterone  receptor-positive,  and  25%  of 
ER+/progesterone  receptor-positive  tumors  (6). 

The  ability  of  IRF-1  to  induce  growth  arrest  is  associated  with  the 
induction  of  various  genes  including  p53-dependent  and  -independent 
events  (26,  27)  and  interactions  that  may  include  p2lwafl/c,P1  (26). 
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Fig.  5.  Doxorubicin  induction  of  IRF-1  mRNA  expression.  A,  representative  RNase 
protection  assay.  36B4,  loading  control.  B,  MCF-7  cells  were  treated  with  1  fx m  doxoru¬ 
bicin.  Data  represent  mean  £  SE  of  three  determinations,  in  which  absorbance  is 
expressed  as  a  ratio  of  IRF-1 :36B4.  *,  P  <  0.01;  Dunnett’s  test. 
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Fig.  6.  Doxorubicin  induces  IRF-l  mRNA  expression  in  both  ER+  and  ER-  breast 
cancer  cell  lines.  RNasc  protection  assay  of  ER+  MCF-7,  T47D,  MCF-7/LCC9,  and 
ZR75  cells  and  ER—  MDA-MB-231  cells  treated  with  1  fx m  doxorubicin  ( DX)  for  24  h. 
A,  representative  data  from  RNasc  protection  assay;  D,  graphical  representation  of  RNase 
protection  assays  of  ER+  and  ER—  breast  cancer  cel!  lines  treated  with  doxorubicin.  Data 
represent  mean  ±  SE  of  three  independent  determinations,  in  which  absorbance  is 
expressed  as  a  ratio  of  IRF-1 :36B4  and  represented  as  a  percentage  of  vehicle  control- 
trcaled  cells.  **,  P  0.001  for  doxorubicin  versus  control-treated  cells  within  each  cell 
line  for  MCF-7  and  MCF-7/LCC9  cells;  *,  P  <  0.05  for  MDA-MB-231  cells;  Student’s 
t  test. 
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Fig.  7.  Inhibition  of  basal  and  IFN-y  stimulated  ISRE  activity  by  dnTRF-1  in  MCF-7 
cells.  Data  represent  mean  ±  SE  (representative  experiment  of  four  independent  repli¬ 
cates),  in  which  data  are  represented  by  relative  light  units.  Where  appropriate,  MCF-7 
cells  were  treated  with  500  IU  of  IFN-y.  *,  P  ^  0.001  for  all  transfections  versus  control; 
Dunnelt’s  test. 

Interactions  requiring  both  p21WI,n/clpl  and  p53  may  not  be  central 
components  in  antiestrogen  signaling  through  IRF-l.  Nonetheless, 
preliminary  data  suggest  an  increase  in  p21v'“n/c,pl  mRNA  expression 
after  ICI  182,780  treatment  in  MCF-7  and  T47D  cells  (2.91  ±  0.89- 
fold),  consistent  with  both  a  previous  report  on  p2 1  wf,ri/c,p|  regulation 
by  ICI  182,780  (56)  and  activation  of  IRF-l.  MCF-7  cells  express 
wild-type  p53,  whereas  T47D  cells  express  a  mutant  and  nonfunc¬ 
tional  p53  (57),  but  both  cell  lines  are  responsive  to  antiestrogen- 
induced  apoptosis  in  a  manner  that  remains  sensitive  to  the  effects  of 
dnIRF-1.  However,  a  role  for  p53/p21wafl/clpl  signaling  in  the  cell 
cycle  effects  of  antiestrogens  cannot  be  excluded. 


The  ability  of  ICI  182,780  to  induce  apoptosis  through  IRF-l 
activity  is  likely  mediated  through  changes  in  caspase  activation. 
IRF-l  can  induce  several  caspases  (27,  29,  30),  and  inhibition  of 
caspase  activity  blocks  antiestrogen-induced  apoptosis  (58).  A  spe¬ 
cific  requirement  for  caspase-3  seems  unlikely  because  this  caspase  is 
not  expressed  in  MCF-7  cells  (59).  IRF-l  signaling  through  caspase-1 
(27),  caspase-7  (29),  and  caspase-8  (30)  is  strongly  implicated.  For 
example,  IRF-l  induces  caspase-1  (60),  which  can  regulate  apoptosis 
in  normal  mammary  epithelial  cells  (61).  Overexpression  of  caspase-1 
is  lethal  in  MCF-7  cells  (62).  Caspase-7  is  expressed  in  MCF-7  cells 
and  may  substitute  for  the  loss  of  caspase-3  in  these  cells  (63).  IFN-y, 
which  induces  IRF-l  activity  in  MCF-7  cells  (Fig.  7),  is  reported  to 
sensitize  both  MCF-7  and  MDA-MB-231  cells  to  apoptosis  through 
inducing  caspase-8  (64).  TAM  can  induce  caspase-8  (58),  and  con¬ 
sistent  with  our  observations  in  T47D  cells,  caspase-8-induced  apo¬ 
ptosis  occurs  independent  of  p53  (64).  Studies  to  determine  which 
caspases  are  functionally  involved  in  IRF-l  signaling  in  breast  cancer 
are  currently  in  progress. 

Although  data  in  this  study  show  reduced  IRF-l  expression  and  loss 
of  its  endocrine  regulation  in  antiestrogen-resistant  cells,  the  level  of 
IRF-l  activity  in  cells  is  also  affected  by  protein-protein  interactions 
with  the  nucleolar  phosphoprotein  nuclcophosmin  (NPM;  Ref.  65). 
Not  previously  reported  in  breast  cancer  cells,  we  now  have  prelim¬ 
inary  data  to  suggest  that  this  functional  interaction  also  occurs  in 
T47D  cells  (not  shown).  NPM  is  an  estrogen-induced  protein  in 
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Fig.  8.  Inhibition  of  ICI  182,780  effects  on  cell  proliferation  by  dnIRF-1,  A.  MCF-7 
cells  with  and  without  constitutive  dnIRF-1  expression.  B,  T47D  cells  with  and  without 
constitutive  dnIRF-i  expression.  For  each  cell  line,  cells  were  treated  with  ethanol  vehicle, 
100  nM  ICI  182,780  or  1  jam  ICI  182,780  for  3  days.  Data  represent  mean  ±  SE  of  three 
determinations.  *,  P  =  0.002,  MCF-7  100  nM  ICI  182,780;  P  =  0.013,  MCF-7  I  jam  ICI 
182,780;  Student’s  /  test.  P  =  0.002,  T47D  100  nM  ICI  182,780;  P  =  0.043,  T47D  1  jam 
ICI  182,780;  Student’s  /  test. 
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MCF-7  cells  that  is  down-regulated  by  antiestrogens  (66),  and  its 
expression  is  increased  in  MCF-7/LCC9  when  compared  with  MCF- 
7/LCC1  cells  (19).  Thus,  in  addition  to  down-regulating  1RF-1  mRNA 
expression,  antiestrogen-resistant  cells  have  up-regulated  expression 
of  an  endogenous  inhibitor  (NPM).  Interestingly,  we  have  previously 
shown  that  NPM  autoantibody  levels  are  lower  in  TAM-trcated  pa¬ 
tients,  suggesting  that  NPM/IRF-1  interactions  also  may  be  clinically 
relevant  (67). 

it  seems  likely  that  an  acquired  antiestrogen  resistance  phenotype  is 
conferred  not  by  the  alteration  of  a  single  gene  or  signal  transduction 
pathway  but  rather  through  the  perturbation  of  a  signaling  network  of 
integrated  signaling  pathways  (15).  The  data  presented  here  are  con¬ 
sistent  with  cell  signaling  through  1RF-1  being  a  key  component  or 
node  in  such  a  signaling  network.  Activity  as  a  signaling  node  is 
implied  by  (a)  the  potential  for  diversity/redundancy  of  signaling  to  a 
key  end  point  (apoptosis)  downstream  of  IRF-1  (cooperation  with 
p53,  p21  w"n/cipl,  and  regulation  of  several  caspases);  (b)  the  redun¬ 
dancy  apparent  in  regulating  IRF-1  activity  (down-regulation  of  basal 
transcription,  loss  of  ER-mediated  transcription,  and  concurrent  up- 
regulation  of  the  endogenous  inhibitor  NPM);  (c)  the  apparent  spec¬ 
ificity  for  antiestrogen  resistance  (cytotoxic  drugs  can  induce  IRF-1  in 
antiestrogen-resistant  cells);  and  ( d)  the  altered  regulation  of  IRF-1 
activity/expression  in  models  of  both  de  novo  and  acquired  antiestro¬ 
gen  resistance.  This  node  may  be  important  in  affecting  other  signals 
in  breast  cancer  cells.  For  example,  IRF-1  is  downstream  of  tumor 


T47D 

Fig.  9.  Expression  of  dnIRF-1  does  not  affect  ICI  182,780-induced  changes  in  cell 
cycle  distribution.  A,  MCF-7  cells  with  and  without  constitutive  dnIRF-1  expression;  B, 
T47D  cells  with  and  without  constitutive  dnIRF-1  expression.  For  each  cell  line,  cells 
were  treated  with  either  ethanol  vehicle  or  1 00  nM  ICI  1 82,780  for  3  days.  Data  represent 
mean  ±  SE  of  three  determinations.  *,  P  =  0.40,  MCF-7  0  nM  ICI  182,780  (vehicle 
control);  P  =  0.237,  MCF-7  100  nM  ICI  182,780;  Student’s  t  test.  P  =  0.962,  T47D  0  nM 
ICI  182,780;  P  =  0.836,  T47D  100  nM;  Student’s  ( test. 
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Fig.  10.  ICI  182,780-induced  apoptosis  is  reduced  by  dnIRF-l  expression.  MCF-7  cells 
with  and  without  constitutive  dnlRF-1  expression.  Data  represent  mean  ±  SE  of  three 
independent  replicate  experiments.  dnlRF-1-  and  control-transfected  cells  were  treated 
with  either  ethanol  vehicle  or  100  nvi  ICI  182,780  for  3  days.  Annexin  V  analysis  was 
measured  by  FACS.  *,  P  <  0.05  for  ICI  1 82,7 80-trea ted  control  transfectants  versus 
dnIRF-1  transfectants  (Dunnett’s  test). 


necrosis  factor  signaling,  and  both  tumor  necrosis  factor  a  and  its 
receptor  tumor  necrosis  factor  R1  are  down-regulated  in  MCF7/LCC9 
cells,  implying  a  cross-resistance  to  tumor  necrosis  factor-mediated 
events  (15,  19). 

The  ability  of  doxorubicin  to  induce  IRF-1  in  antiestrogen-resistant 
cells  and  to  inhibit  proliferation  in  these  cells  is  a  clinically  relevant 
phenotype.  We  and  others  (36,  68)  have  detected  IRF-1  expression  by 
immunohistochcmistiy  in  breast  cancer  specimens,  this  pattern  of 
expression  being  consistent  with  a  potential  tumor  suppressor  role  for 
IRF-1  (2 1 , 36).  Our  preliminary  data  suggest  that  the  pattern  of  IRF-1 
expression  in  breast  cancers,  as  measured  by  immunohistochemistry, 
is  consistent  with  other  components  of  our  network.  For  example,  we 
detect  an  inverse  pattern  of  expression  between  IRF-1  and  nuclear 
factor  kB  as  seen  in  MCF7/LCC1  versus  MC.F7/LCC9  cells  (19,  68). 
These  observations  may  ultimately  lead  to  a  better  ability  to  identify 
patients  that  will  respond  to  antiestrogens  and  to  predict  which  pa¬ 
tients  will  ultimately  develop  antiestrogen  resistance.  Interfering  with 
the  putative  “IRF-1  node”  may  allow  for  the  development  of  novel 
therapeutic  strategies  in  endocrine-resistant  breast  cancers. 
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The  nuclear  factor  kB  inhibitor  parthenolide  restores 
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Abstract 

The  molecular  mechanisms  underlying  the  acquisition  of 
resistance  to  the  antiestrogen  Faslodex  are  poorly 
understood,  although  enhanced  expression  and  activity 
of  nuclear  factor  k,B  (NFkB)  have  been  implicated  as  a 
critical  element  of  this  phenotype.  The  purpose  of  this 
study  was  to  elucidate  the  mechanism  by  which  NFkB 
up-regulation  contributes  to  Faslodex  resistance  and  to 
determine  whether  pharmacologic  inhibition  of  NFkB  by 
the  small  molecule  parthenolide  could  restore  Faslodex- 
mediated  suppression  of  cell  growth.  Basal  expression  of 
multiple  NFKB-related  molecules  in  MCF7-derived  LCC1 
(antiestrogen-sensitive)  and  LCC9  (antiestrogen-resistant) 
breast  cancer  cells  was  determined,  and  cells  were  treated 
with  Faslodex  or  parthenolide.  The  effect  of  these  drugs 
either  singly  or  in  combination  was  assessed  by  cell 
proliferation,  estrogen  receptor  (ER) -dependent  tran¬ 
scriptional  activation,  cell  cycle  analysis,  and  apoptosis 
assays.  Expression  of  the  p65  NFkB  subunit  and  the 
upstream  NFkB  regulator  IkB  kinase  7/NFkB  essen¬ 
tial  modulator  were  increased  in  the  resistant  MCF7/ 
LCC9  cells  (P  =  0.001  and  0.04,  respectively).  Whereas 
MCF7/LCC9  cells  were  unresponsive  to  Faslodex  alone, 
parthenolide  effectively  inhibited  MCF7/LCC9  cell  prolif¬ 
eration  and  the  combination  of  Faslodex  and  parthenolide 
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resulted  in  a  4-fold  synergistic  reduction  in  cell  growth 
(P  =  0.03).  This  corresponded  to  a  restoration  of 
Faslodex-induced  apoptosis  (P  =  0.001),  with  no  observ¬ 
able  changes  in  ER-dependent  transcription  or  cell  cycle 
phase  distribution.  Because  parthenolide  has  shown 
safety  in  Phase  I  clinical  trials,  these  findings  have  direct 
clinical  relevance  and  provide  support  for  the  design  of 
clinical  studies  combining  antiestrogens  and  parthenolide 
in  ER-positive  breast  cancer.  [Mol  Cancer  Ther  2005; 
4(11:33-41] 


Introduction 

Antiestrogens  inhibit  the  function  of  the  estrogen  receptor 
(ER),  a  nuclear  transcription  factor  that  directs  the 
expression  of  genes  that  contribute  to  proliferation  and 
cell  growth  (reviewed  in  refs.  1,  2).  The  most  frequently 
prescribed  is  the  nonsteroidal  antiestrogen  tamoxifen, 
which  has  been  shown  to  be  highly  effective  in  both  the 
treatment  of  ER-positive  breast  tumors  and  in  reducing 
breast  cancer  incidence  in  women  at  high  risk  for  the 
disease.  However,  most  ER-positive  tumors  become  estro¬ 
gen  independent  and  develop  resistance  to  antiestrogen 
therapy,  whereas  the  remainder  ( —  30%)  exhibit  de  novo  or 
intrinsic  resistance.  Once  resistance  has  developed,  treat¬ 
ment  with  most  nonsteroidal  antiestrogens  is  usually 
unsuccessful. 

In  contrast,  the  steroidal  antiestrogen  Faslodex  (ICI 
182,780;  ICI)  induces  significant  clinical  responses  in 
patients  whose  tumors  have  acquired  tamoxifen  resistance 
(3).  The  effectiveness  of  Faslodex  in  patients  with  tamox¬ 
ifen-resistant  disease  is  similar  to  that  of  the  aromatase 
inhibitor  anastrozole,  and  several  clinical  trials  have  shown 
that  Faslodex  may  be  a  viable  alternative  to  nonsteroidal 
antiestrogens  and  aromatase  inhibitors  as  a  first-line 
endocrine  treatment  (4).  Faslodex  stimulates  degradation 
of  the  ER  and  prevents  receptor  dimerization,  inhibiting 
estrogen-dependent  gene  transcription  (5,  6).  As  a  pure 
antagonist  of  the  ER,  Faslodex  is  not  associated  with  the 
increased  risk  for  endometrial  cancer  that  is  seen  with 
tamoxifen  (7). 

The  antiestrogen  resistance  phenotype  is  complex,  in¬ 
volving  many  changes  at  the  cellular  and  molecular  levels. 
Antiestrogens  are  cytostatic,  inducing  a  Go-G}  block  in 
breast  cancer  cells  in  culture  (1,8).  However,  these  drugs  are 
also  capable  of  actively  inducing  programmed  cell  death  or 
apoptosis,  which  is  consistent  with  the  ability  of  anties¬ 
trogens  to  increase  overall  survival  (9).  One  way  in  which 
breast  cancer  cells  may  become  antiestrogen  resistant  is 
through  changes  in  gene  networks  that  control  cell 
proliferation  and  apoptosis  (10).  To  test  this  hypothesis, 
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we  developed  several  variant  cell  lines  from  the  estrogen- 
dependent  and  antiestrogen-sensitive  MCF-7  breast  cancer 
cells  (11,  12).  MCF7/LCC1  cells  are  estrogen  independent 
but  remain  responsive  to  antiestrogens;  MCF7/LCC9  cells 
are  derivatives  of  MCF7/LCC1  that  have  acquired  resis¬ 
tance  to  Faslodex.  Similar  to  what  has  been  observed  in 
breast  cancer  patients,  MCF7/LCC9  cells  are  cross-resistant 
to  the  nonsteroidal  antiestrogen  tamoxifen  (2). 

Several  genes  were  found  to  be  altered  in  the  resistant 
MCF7/LCC9  cells,  when  their  transcriptomes  were  com¬ 
pared  with  that  of  their  antiestrogen-sensitive  MCF7/LCC1 
parental  cells  by  serial  analysis  of  gene  expression  and 
microarray  analysis  (10).  For  example,  we  implicated  loss 
of  the  putative  tumor  suppressor  interferon  regulatory 
factor-1  (IKFI)  in  acquired  resistance  and  have  recently 
shown  IRFI  to  be  a  key  mediator  of  the  proapoptotic  effects 
of  Faslodex  in  MCF-7  cells  (13). 

Altered  expression  of  the  p65/RelA  member  of  the 
nuclear  factor  kB  (NFkB)  transcription  factor  family,  which 
can  form  functional  heterodimers  with  IFN  regulatory 
factor-1  (14),  also  was  strongly  implicated  in  acquired 
Faslodex  resistance.  mRNA  levels  of  p65/RelA  are  up- 
regulated  2-fold  in  the  MCF7/LCC9  cells,  NFKB-depcndent 
transcription  are  increased  10-fold,  and  MCF7/LCC9  cells 
exhibit  a  greater  sensitivity  to  the  growth  inhibitory  effects 
of  parthenolide,  a  small  molecule  inhibitor  of  NFkB  (10). 
These  data  strongly  but  indirectly  implicate  NFkB  action  in 
acquired  antiestrogen  resistance. 

The  NFkB  family  contains  five  members  that  form 
dimers  and  regulate  the  transcription  of  various  genes 
including  cytokines,  cell  adhesion  molecules,  the  pro- 
proliferative  proteins  c-myc  and  cyclin  Dl,  and  several 
inhibitors  of  apoptosis  (15).  Inhibitors  of  the  NFkB  pathway 
show  promise  as  anticancer  and  anti-inflammatory  agents 
(16).  Parthenolide,  a  sesquiterpene  lactone  that  was  first 
isolated  from  the  feverfew  herb  (Tamcetum  pnrthenium ) 
native  to  Central  America  (17),  is  a  relatively  specific  small 
molecule  inhibitor  of  NFkB  (18).  Parthenolide  and  other 
members  of  the  sesquiterpene  lactone  class  have  garnered 
recent  attention  as  promising  candidates  for  cancer 
treatment  either  as  single  agents  or  in  combination  with 
other  cytotoxic  drugs  (19,  20).  For  example,  parthenolide 
has  anti-inflammatory,  anticancer,  and  antiangiogenic 
properties  and  has  successfully  undergone  phase  I/II 
clinical  trials  (21,  22). 

Constitutive  NFkB  activity  is  widely  observed  in  many 
tumor  types  (23),  including  breast  cancer  where  it  is 
associated  with  resistance  to  apoptosis-inducing  agents 
(24).  In  many  tumor  lines,  autocrine  secretion  of  cytokines 
and  growth  factors  has  recently  been  implicated  in  the 
constitutive  activation  of  NFkB  (25).  Importantly,  NFkB 
activity  also  increases  in  breast  cancer  cells  as  they  acquire 
the  ability  to  grow  in  the  absence  of  estrogen  (26,  27).  These 
findings  strongly  implicate  NFkB  signaling  in  the  control  of 
breast  cancer  cell  growth  and  response  to  antiestrogens. 

In  this  study,  we  sought  to  clarify  the  mechanism  by 
which  NFkB  up-regulation  may  affect  resistance  to  Faslo¬ 
dex  and  determine  whether  pharmacologic  inhibition  of 


NFkB  could  restore  sensitivity  to  the  drug.  We  show  here 
that  in  addition  to  p65/RelA,  expression  of  the  upstream 
regulator  NFkB  essential  modulator/kB  kinase  -y  (NEMO/ 
IKKy)  is  also  increased  in  the  resistant  cells.  The  NFkB 
inhibitor  parthenolide  efficiently  inhibits  cell  growth  and 
restores  sensitivity  to  Faslodex  by  synergistically  enhancing 
apoptosis.  Our  data  indicate  that  inhibition  of  NFkB  may  be 
a  successful  approach  in  the  treatment  of  ER-positive  breast 
cancers  that  have  acquired  resistance  to  antiestrogen 
therapy.  NFkB  inhibition  also  may  reduce  the  incidence 
or  delay  the  onset  of  antiestrogen  resistance.  These  data 
provide  support  for  considering  the  design  of  clinical 
studies  combining  antiestrogens  and  parthenolide  in  ER+ 
breast  cancer. 


Materials  and  Methods 

Cell  Culture  and  Reagents 

MCF-7 -derived  MCF7/LCC1  and  MCF7/LCC9  cells 
(11,  12)  were  routinely  cultured  in  phenol  red-free  im¬ 
proved  minimal  essential  media  (IMEM;  Biofluids,  Rock¬ 
ville,  MD)  supplemented  with  5%  charcoal-stripped  calf 
serum  (CCS;  CCS-IMEM).  Cells  were  maintained  in  a 
humidified  atmosphere  at  37°C  and  95%  air/5%  CO2.  17fV- 
Estradiol  (estradiol,  E2)  and  parthenolide  were  purchased 
from  Sigma  (St.  Louis,  MO),  and  ICI  182,780  (IC1,  Faslodex) 
was  a  kind  gift  of  Dr.  Alan  Wakeling  (AstraZeneca, 
Macclesfield,  Cheshire,  United  Kingdom). 

Cell  Lysis,  Immunoblotting,  and  Immunoprecipitation 

Cells  were  grown  in  either  10-cm2  dishes  or  T-75  cm2 
tissue  culture  flasks  before  lysis.  To  determine  the  effects  of 
parthenolide  and  1C1  182,780  on  protein  expression,  cells 
were  treated  with  vehicle,  100  nmol/L  ICI  182,780  (the  IC5() 
for  the  control/parental  LCC1  cells),  or  600  nmol/L 
parthenolide  singly  or  in  combination  in  CCS-IMEM  for 
72  hours.  Cells  were  then  lysed  in  modified  radioimmuno- 
precipitation  assay  buffer  [150  mmol/L  NaCl,  50  mmol/L 
Tris  (pH  7.5),  1%  Igepal  CA-630,  and  0.5%  deoxycholate] 
supplemented  with  Complete  Mini  protease  inhibitor 
cocktail  tablets  (Roche,  Mannheim,  Germany)  and 
1  mmol/L  sodium  orthovanadate  phosphatase  inhibitor 
(Sigma).  Lysates  were  clarified  by  centrifugation  and  total 
protein  was  quantitated  using  the  bicinchoninic  acid  assay 
purchased  from  Pierce  (Rockford,  IL).  Whole  cell  lysate 
(20  pg)  was  resolved  by  PAGE  using  NuPAGE  12%  precast 
gels  (Invitrogen,  Carlsbad,  CA).  Proteins  were  then 
transferred  to  nitrocellulose  membranes,  which  were 
probed  w'ith  the  following  antibodies  overnight  at  4°C: 
p65  NFkB  sc-109  (1:800;  Santa  Cruz  Biotechnology,  Santa 
Cruz,  CA),  p50  NFkB  sc-8414  (1:200;  Santa  Cruz  Biotech¬ 
nology),  p52  NFkB  (1:200;  Upstate  Biotechnology,  Char¬ 
lottesville,  VA),  IKKy/NEMO  sc-8330  (1:200;  Santa  Cruz 
Biotechnology),  kBa  sc-371  (1:200;  Santa  Cruz  Biotechnol¬ 
ogy),  phospho-Akt  (Ser473;  1:1,000,  Cell  Signaling,  Beverly, 
MA),  or  Akt  (1:1,000,  Cell  Signaling).  Membranes  were  then 
incubated  with  horseradish  peroxidase-conjugated  sec¬ 
ondary  antibodies  (Amersham  Biosciences,  Piscataway,  NJ) 
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for  1  hour  at  room  temperature  before  enhanced  chemilu¬ 
minescence  (Amersham  Biosciences)  and  exposure  to  film. 
To  confirm  equal  loading  of  the  gels,  membranes  were 
reprobed  with  antibodies  for  p-actin  (1:5,000;  Sigma). 

For  immunoprecipitations,  400  pg  of  cell  lysate  were 
incubated  with  2.5  gL  of  p65  NFkB  antibodies  overnight  at 
4°C  with  rotation.  The  following  day,  30  gL  of  Protein  A- 
Sepharose  beads  (Amersham  Biosciences)  were  added  for  1 
hour  at  4°C  to  recover  the  immune  complexes,  which  were 
then  washed  twice  in  modified  radioimmunoprecipitation 
assay  buffer,  twice  in  Tris-saline  [TN;  50  mmol/L  Tris  (pH 
7.5)  and  150  mmol/L  NaCl],  and  resuspended  in2x  Laemmli 
sample  buffer  before  electrophoresis  as  described  above. 

Cell  Proliferation  Assays 

MCF7/LCC1  and  MCF7/LCC9  cells  were  seeded  at  a 
density  of  1  to  2  x  104  cells  per  well  in  24-well  plates,  and 
24  hours  later  were  treated  with  the  indicated  concentrations 
of  drug  in  CCS-IMEM.  Cells  were  incubated  with  the  drugs 
for  7  days,  and  the  media  were  changed  on  days  3  and  5. 
Cells  were  then  trypsinized,  resuspended  in  PBS  (Biofluids), 
and  counted  using  a  Z1  Single  Coulter  Counter  (Beckman 
Coulter,  Miami,  FL).  At  least  three  independent  experiments 
were  done  in  quadruplicate,  and  data  were  normalized  to 
vehicle-treated  cells.  Data  are  presented  as  the  mean  ±  SE  for 
a  representative  experiment. 

Transcriptional  Reporter  Assays 

The  estrogen  response  element-containing  3xERE-tk-luc 
reporter  plasmid  was  purchased  from  Promega  (Madison, 
WI).  MCF7/LCC1  and  MCF7/LCC9  cells  were  seeded  into 
12-well  plates  at  a  density  of  7  to  8  x  104  cells  per  well.  The 
following  day,  cells  were  transfected  with  0.4  gg  of  luciferase 
reporter  plasmid  and  0.1  gg  pCMV-Renilln  (Promega)  per 
well  using  the  FuGENE  6  transfection  reagent  (Roche, 
Indianapolis,  IN).  Three  hours  post-transfection,  media 
were  changed  and  cells  were  treated  with  100  nmol/L  ICI 
182,780  and/or  600  nmol/L  parthenolide  in  CCS-IMEM  for 
24  hours.  Subsequently,  cells  were  lysed  and  activation  of 
the  luciferase  constructs  was  measured  using  the  Dual 
Luciferase  Assay  Kit  (Promega).  Luminescence  was  quan¬ 
tified  using  a  Lumat  LB  9501  luminometer  (EG&G  Berthold, 
Bundoora  VIC,  Australia).  Luciferase  values  were  normal¬ 
ized  to  Renilln  luminescence,  and  four  independent  experi¬ 
ments  were  done  each  at  least  in  quadruplicate.  Data  are 
presented  as  the  mean  ±  SE  for  all  experiments. 

Cell  Cycle  Assays 

Cells  (n  =  5  x  105)  were  seeded  into  10-cm2  dishes  1  day 
before  treatment  with  100  nmol/L  ICI  182,780  and/or  600 
nmol/L  parthenolide  in  CCS-IMEM  for  24  hours.  Cells 
were  then  analyzed  for  alterations  in  cell  cycle  via 
fluorescence  activated  cell  sorting,  which  was  done  by  the 
Lombardi  Comprehensive  Cancer  Center  Flow  Cytometry 
Shared  Resource  according  to  the  method  of  Vindelov  et  al. 
(28).  Data  are  presented  as  the  mean  ±  SE  for  three 
independent  experiments. 

Apoptosis  Assays 

Cells  (n  =  1  x  105)  were  seeded  onto  18  x  18  mm  glass 
coverslips  in  each  well  of  a  6-well  plate  in  duplicate  and 
the  following  day  were  treated  with  100  nmol/L  ICI 


182,780  and/or  600  nmol/L  parthenolide  in  CCS-IMEM 
for  24  hours.  Cells  were  then  fixed  with  3.7%  formalin  in 
PBS  for  20  minutes  at  room  temperature  prior  to  Annexin 
V  and  propidium  iodide  staining  using  the  Vybrant 
Apoptosis  Assay  Kit  3  purchased  from  Vector  Laborato¬ 
ries  (Burlingame,  CA).  Coverslips  were  then  mounted  on 
glass  slides  using  VectaShield  fluorescence  mounting 
medium  (Vector  Laboratories).  Cells  were  visualized  on 
a  Nikon  E600  fluorescence  microscope  (provided  by  the 
Lombardi  Comprehensive  Cancer  Center  Microscopy 
Shared  Resource),  and  several  random  fields  (3;200  cells) 
were  scored  per  treatment  condition.  The  number  of  cells 
stained  red  (propidium  iodide,  indicating  necrosis)  was 
subtracted  from  the  number  of  cells  stained  green 
(Annexin  V-FITC,  indicating  apoptosis),  and  subsequently 
divided  by  the  total  number  of  cells  seen  by  phase- 
contrast.  Data  are  presented  as  the  percentage  of  apoptotie 
cells  and  represent  the  mean  ±  SE  for  three  independent 
experiments. 

Statistical  Analyses 

Two-tailed  Student's  f  tests  were  used  for  the  compar¬ 
ison  of  two  groups  for  immunoblot,  cell  proliferation,  and 
apoptosis  assays  as  indicated.  For  luciferase  reporter 
assays,  Dunn's  post  hoc  t  test  was  used  to  compare  all 
treatment  groups  following  one-way  ANOVA.  Defining  the 
nature  of  the  interaction  between  Faslodex  and  partheno¬ 
lide  was  done  by  determining  the  R  index  (RI;  ref.  29).  RI 
values  were  obtained  by  calculating  the  expected  cell 
survival  (Scxp;  the  product  of  survival  obtained  with  drug 
A  alone  and  the  survival  obtained  with  drug  B  alone)  and 
dividing  Sexp  by  the  observed  cell  survival  in  the  presence 
of  both  drugs  (Sobs).  Sexp/Sobs  >  1.0  indicates  a  synergistic 
interaction.  This  method  is  an  appropriate  way  to  define 
synergy  in  this  case  because  clinically  relevant  concen¬ 
trations  of  Faslodex  are  ineffective  on  cell  proliferation  or 
apoptosis  in  the  resistant  MCF7/LCC9  cells  when  given  as 
a  single  agent  (29). 

Results 

p65/RelA  and  NEMO/IKKy  Are  Up-Regulated  in 
Antiestrogen-Resistant  Cells 

Our  previous  studies  identified  a  2-fold  up-regulation  of 
p65/RelA  mRNA  in  Faslodex-resistant  MCF7/LCC9  cells 
by  expression  microarray  analysis  (10).  To  confirm  altered 
expression  of  p65  at  the  protein  level  and  to  examine  other 
NFkB  family  members  and  regulatory  molecules,  whole 
cell  lysates  were  prepared  from  MCF7/LCC9  cells  and 
antiestrogen-sensitive  MCF7/LCC1  cells  and  subjected  to 
SDS-PAGE  and  immunoblot  analysis.  Similar  to  mRNA 
levels,  p65/RelA  protein  is  increased  ~  2-fold  in  the 
MCF7/LCC9  cells  (Fig.  1A,  P  =  0.001).  In  contrast,  expres¬ 
sion  of  the  p50  subunit  of  the  NFkB  heterodimer  (Fig.  IB, 
P  =  0.35)  or  of  p52  NFkB2  (data  not  shown)  is  not  different 
between  the  cell  lines. 

Transcriptional  activity  of  the  p65/p50  heterodimer  is 
modulated  by  the  inhibitor  IkB,  which  is  in  turn  negatively 
regulated  by  the  IKK  complex  comprised  of  IKKa,  IKKp, 
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and  the  scaffolding  protein  NEMO/IKKy  (30).  To  deter¬ 
mine  whether  MCF7/LCC9  cells  exhibit  changes  in  these 
regulatory  molecules,  lysates  were  immunoblotted  for 
NEMO/IKKy  and  IkBct  (Fig.  1C  and  D).  Whereas  there  is 
no  significant  change  in  IkBc*  expression  (P  =  0.10),  a 
significant  2-fold  increase  in  the  level  of  NEMO/IKK-y  is 
apparent  in  MCF7/LCC9  cells  (P  =  0.04).  NEMO/IKKy  is 
required  for  activity  of  the  IKK  complex  and  the  inhibitory 
phosphorylation  of  IkB  in  response  to  inflammatory  stimuli 
that  activate  NFkB  (31),  and  dysregulation  of  NEMO  is 
linked  to  several  human  pathologies  (32).  These  data 
suggest  that  NEMO  may  also  play  a  role  in  the  response 
of  breast  cancer  to  antiestrogens. 

To  examine  whether  the  binding  of  p65  and  p50  was 
altered  in  antiestrogen-resistant  cells,  cell  lysates  were 
immunoprecipitated  with  p65  antibodies  and  immune 
complexes  were  captured  and  subjected  to  SDS-PAGE  as 
described  above  (Fig.  IE).  No  clear  differences  in  p65/p50 
complex  formation  were  found  between  MCF7/LCC1  and 
MCF7/LCC9  cells. 

Independent  of  the  IKK-IkB  signaling  pathway,  NFkB 
can  also  be  activated  by  phosphatidylinositol  3-kinase 
(PI3K);  PI3K-mediated  activation  of  Akt  can  enhance 
NFkB  transcriptional  activity  without  the  degradation  of 
IkB  (33).  Because  overexpression  of  active  Akt  has  also 
been  shown  to  induce  resistance  to  antiestrogens  and 
cytotoxic  drugs  (34),  MCF7/LCC1  and  MCF7/LCC9  cell 
lysates  described  above  were  immunoblotted  for  phospho- 
Serine  473  Akt  (Fig.  IF).  No  difference  in  the  level  of 
activated  phospho-Akt  is  observed  in  the  antiestrogen- 
resistant  MCF7/LCC9  cells,  suggesting  that  Akt  — *  NFkB 
signaling  is  not  the  only  pathway  through  which  cells  can 
modulate  NFkB  activation  and  acquire  resistance  to 
Faslodex. 

Inhibition  of  NFkB  by  Parthenolide  Restores  Faslodex 
Sensitivity  to  MCF7/LCC9  Cells 

We  have  previously  reported  that  MCF7/LCC9  cells 
are  more  sensitive  than  MCF7/LCC1  cells  to  growth 
inhibition  by  parthenolide,  suggesting  that  these  cells,  in 
which  p65/RelA  is  up-regulated,  are  more  dependent  on 
NFkB-cI riven  cell  growth  (10).  Faslodex  (100  nrnol/L) 
approximates  the  ICso  for  proliferation  in  antiestrogen- 
sensitive  MCF7/LCC1  cells  but  is  ineffective  in  MCF7/ 
LCC9  cells  (Fig.  2A,  P  =  0.01).  To  determine  whether 


Figure  1 .  Expression  of  NFkB  family  members  and  upstream  regulatory 
molecules.  A-D,  quantitiation  and  representative  immunoblots  of  p65/ 
RelA,  p50  NFkB,  NEMO/IKK7,  and  IkBu  levels  in  MCF7/LCC1  and  MCF7/ 
LCC9  cells.  Lysates  (20  pg)  were  separated  by  SDS-PAGE,  transferred  to 
nitrocellulose  membrane,  and  immunoblotted.  /i-actin,  loading  control. 
Columns,  mean  of  at  least  3  independent  experiments;  bar,  ±SE.  Ps  were 
calculated  by  Student's  t  test.  E,  coimmunoprecipitation  of  p65  and  p50. 
Lysates  (400  pg)  were  immunoprecipitated  with  polyclonal  anti-p65 
antibodies;  immune  complexes  were  isolated  and  separated  by  SDS- 
PAGE,  transferred  to  nitrocellulose,  and  immunoblotted.  F,  Akt  activity  is 
not  altered  in  antiestrogen-resistant  MCF7/LCC9  cells.  Lysates  (20  pg) 
were  separated  by  SDS-PAGE,  transferred  to  nitrocellulose,  and  immuno¬ 
blotted  with  antibodies  specific  for  phospho-Ser473  of  Akt.  The  mem¬ 
brane  was  then  stripped  and  reprobed  for  total  Akt. 
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Figure  2.  Parthenolide  inhibits  the  proliferation  of  antiestrogen-resistant 
cells  and  partially  restores  Faslodex  sensitivity.  A,  MCF7/LCC9  cells  are 
unresponsive  to  Faslodex.  Cells  were  seeded  in  quadruplicate  and  treated 
with  1 00  nmol/L  Faslodex  in  CCS-IMEM  for  7  d  before  counting.  Columns, 
mean  from  a  single  representative  experiment  of  relative  proliferation 
(relative  to  vehicle-treated  control);  bar,  ±SE.  P  was  calculated  by 
Student's  t  test.  The  experiment  was  independently  done  at  least  thrice. 
B,  Faslodex  and  parthenolide  synergistically  inhibit  MCF7/LCC9  cell 
proliferation.  Cells  were  seeded  in  quadruplicate  and  treated  with  0  to 
6,000  nmol/L  parthenolide  in  the  presence  or  absence  of  100  nmol/L 
Faslodex  in  CCS-IMEM  for  6  d.  Points,  mean  of  relative  proliferation;  bars, 
±SE.  P  -  0.034  versus  100  nmol/L  parthenolide  without  Faslodex  by 
Student's  f  test;  Rl  =  1.82.  it,  P  -  0.05  versus  600  nmol/L  parthenolide 
without  Faslodex;  Rl  =  1.48. 


inhibition  of  NFkB  activity  could  restore  Faslodex 
sensitivity,  MCF7/LCC9  cells  were  treated  with  increasing 
concentrations  of  parthenolide  in  the  presence  or  absence  of 
100  nmol/L  Faslodex.  In  the  absence  of  Faslodex,  parthe¬ 
nolide  effectively  inhibits  MCF7/LCC9  cell  proliferation 
with  an  IC50  of  500  to  600  nmol/L.  However,  the  addition 
of  Faslodex  generates  a  significant  nearly  5-fold  sensitiza¬ 
tion,  where  50%  growth  inhibition  occurs  at  a  concentration 
of  100  nmol/L  parthenolide  (Fig.  2B;  P  =  0.034  for 
parthenolide  plus  Faslodex  compared  with  parthenolide 
alone). 

The  interaction  of  Faslodex  and  parthenolide  is 
synergistic  in  MCF7/LCC9  cells,  generating  an  Rl  value 
of  1.82.  Treatment  with  100  nmol/L  Faslodex  and  600 
nmol/L  parthenolide  also  produces  a  greater  than 


additive  inhibition  of  cell  proliferation  (P  =  0.05,  Rl  = 
1.48).  These  data  strongly  suggest  that  the  up-regulated 
NFkB  activity  present  in  MCF7/LCC9  cells  is  a  major 
contributor  to  the  antiestrogen  resistance  phenotype. 

Parthenolide  and  Faslodex  Synergistically  Increase 
Apoptosis 

We  subsequently  sought  to  define  the  mechanism  by 
which  parthenolide  and  Faslodex  synergistically  inhibit 
the  growth  of  MCF7/LCC9  cells.  A  primary  action  of 
antiestrogens  is  to  antagonize  endogenous  estrogen  and 
block  ER  function;  Faslodex  can  achieve  this  by  affecting 
receptor  turnover  (5).  We  asked  whether  parthenolide  can 
restore  Faslodex-mediated  inhibition  of  ER-dependent 
transcriptional  activity  (Fig.  3).  MCF7/LCC1  and  MCF7/ 
LCC9  cells  were  cotransfected  with  an  ERE-tk-luciferase 
reporter  vector  and  the  pCMV-Rcnilla  control  vector.  Three 
hours  post-transfection,  cells  were  treated  with  estradiol, 
Faslodex,  and/or  parthenolide  for  24  hours  before 
performing  dual-luciferase  promoter-reporter  assay's. 

MCF7/LCC1  cells  exhibit  a  basal  ERE-luciferase  activity 
that  is  enhanced  8-fold  by  estradiol  treatment  and  almost 
abolished  by  Faslodex.  In  contrast,  MCF7/LCC9  cells 
express  a  higher  basal  ERE-luciferase  activity  that  is 
slightly  enhanced  by  estradiol  but  is  not  inhibited  by 
Faslodex  treatment.  Whereas  transcription  from  an  NFkB- 
dependent  reporter  is  inhibited  by  600  nmol/L  partheno¬ 
lide  in  MCF7/LCC9  cells  (data  not  shown),  parthenolide 
either  alone  or  in  combination  with  Faslodex  has  no 
statistically  significant  effect  on  ERE-luciferase  activity  in 
MCF7/LCC9  cells,  suggesting  that  the  mechanism  of  their 
antiproliferative  synergy  does  hot  involve  tire  regulation 
of  ER-dependent  transcriptional  events. 

Treatment  with  antiestrogens  such  as  Faslodex  can  have 
a  cytostatic  effect  on  cell  growth,  typically  manifested  as 


Figure  3.  Combined  treatment  with  Faslodex  and  parthenolide  does  not 
inhibit  ER-dependent  transcriptional  activity.  MCF7/LCC 1  and  MCF7/LCC9 
cells  were  transfected  in  quadruplicate  with  ERE-tk-luciferase  and  pCMV- 
Renilla  constructs  prior  to  treatment  with  1 0  nmol/L  estradiol,  1 00  nmol/L 
Faslodex,  and  600  nmol/L  parthenolide  singly  or  in  combination  (or  ethanol 
vehicle)  in  CCS-IMEM  for  24  h.  Columns,  mean  of  the  ratio  of  luciferase- 
to-Renilla  activity  (relative  light  units)  for  four  independent  experiments; 
bars,  ±SE.  P  <  0.001  for  all  treatment  groups  by  one-way  ANOVA, 
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Figure  4.  Combined  treatment  with  Faslodex  and  parthenolide  has  no 
effect  on  the  cell  cycle  profile  of  MCF7/LCC9  cells.  Cells  were  treated 
with  100  nmol/L  Faslodex,  600  nmol/L  parthenolide,  Faslodex  + 
parthenolide,  or  ethanol  vehicle  in  CCS-IMEM  for  24  h  before  cell  cycle 
analysis.  Columns,  mean  for  three  independent  experiments  {%  total 
cells);  bars,  ±SE. 

an  accumulation  of  cells  in  the  Go-Gi  phase  of  the  cell 
cycle  (1,  8).  hi  some  cell  systems,  parthenolide  can  arrest 
cells  at  the  G2-M  phase  transition  (35).  To  test  whether 
parthenolide  restored  the  cytostatic  activities  of  Faslodex 
or  induced  a  GrM  blockade,  MCF7/LCC9  cells  were 
treated  with  Faslodex  ±  parthenolide  or  ethanol  vehicle  for 
24  hours  prior  to  cell  cycle  analysis  (Fig.  4).  Parthenolide 
alone  or  in  combination  with  Faslodex  does  not  alter 
the  MCF7/LCC9  cell  cycle  profile,  indicating  that  a  block 
in  cell  cycle  progression  does  not  explain  the  synergistic 
reduction  in  cell  growth. 

Faslodex  and  other  antiestrogens  actively  promote 
apoptosis,  and  parthenolide  has  been  shown  to  cooper¬ 
atively  enhance  apoptosis  induced  by  other  cytotoxic 
agents  such  as  paclitaxel  and  4-hydroxyphenylretinamide 
(19,  20).  Therefore,  we  measured  the  effects  of  Faslodex  ± 
parthenolide  or  ethanol  vehicle  on  apoptosis  as  detected 
by  immunostaining  for  FITC-conjugated  Annexin  V  and 
propidium  iodide  staining  (Table  1).  Approximately  3%  of 
vehicle-treated  and  4%  of  Faslodex-treated  MCF7/LCC9 
cells  undergo  apoptosis.  In  contrast,  parthenolide  treatment 
increases  the  apoptotic  fraction  to  nearly  10%;  upon 
cotreatment  with  Faslodex  and  parthenolide,  18%  of  the 
cells  undergo  apoptosis.  Importantly,  the  level  of  apoptosis 
seen  in  the  presence  of  the  Faslodex/parthenolide  combi¬ 
nation  was  essentially  identical  to  that  induced  by  Faslodex 
alone  in  the  antiestrogen-sensitive  LCC1  cells  (Table  1).  The 
strong  induction  of  apoptosis  in  MCF7/LCC9  cells  seen  in 
the  presence  of  both  drugs  is  statistically  significant  com¬ 
pared  with  either  Faslodex  or  parthenolide  alone  (P  =  0.001 
and  P  =  0.01,  respectively).  The  calculated  RI  =  2.28  for 
the  parthenolide/Faslodex  interaction  indicates  synergistic 
induction  of  apoptosis  in  antiestrogen-resistant  MCF7/ 
LCC9  cells. 

Parthenolide  stabilizes  the  inhibitor  IkB,  leading  to 
the  retention  of  p65  in  the  cytoplasm  in  an  inactive  state 
(36).  Therefore,  we  measured  the  effects  of  Faslodex  ± 


parthenolide  or  ethanol  vehicle  on  kBa  expression  (Fig.  5). 
Since  protein  levels  of  IkBc*  were  unchanged  in  MCF7/ 
LCC9  cells  regardless  of  treatment,  parthenolide  may  be 
acting  through  other  alternative  mechanisms  to  synergize 
with  Faslodex  and  restore  the  apoptotic  response  to 
antiestrogen-resistant  MCF7/LCC9  cells. 

Discussion 

Our  previous  studies  reported  the  p65/RelA  subunit  of 
NFkB  as  being  up-regulated  in  MCF-7-derived  MCF7/ 
LCC9  breast  cancer  cells  that  had  acquired  resistance  to 
Faslodex  (10,  11).  We  have  now  identified  additional 
changes  in  the  expression  of  NFkB  pathway  members  in 
these  cells  and  showed  that  pharmacologic  inhibition  of 
NFkB  restores  Faslodex  sensitivity  by  markedly  enhancing 
apoptosis.  Because  the  NFkB  inhibitor  parthenolide  is 
currently  being  investigated  in  clinical  trials  (21),  these 
findings  have  direct  clinical  relevance  and  provide  support 
for  the  design  of  clinical  studies  combining  antiestrogens 
and  NFkB  inhibitors  such  as  parthenolide  in  ER+  breast 
cancer. 

Protein  expression  of  the  p65/RelA  subunit  of  NFkB  is 
increased  ~  2-fold  in  MCF7/LCC9  cells  when  compared 
with  antiestrogen-sensitive  MCF7/LCC1  cells;  this  agrees 
with  the  up-regulation  in  mRNA  levels  previously  ob¬ 
served  (10).  Flowever,  NFKB-dependent  transcriptional 
activity  is  elevated  almost  10-fold  in  MCF7/LCC9,  imply¬ 
ing  that  other  elements  of  the  NFkB  signaling  pathways  are 
activated  in  these  cells.  We  found  no  changes  in  p50 
expression  or  association  with  p65;  there  were  also  no 
alterations  in  expression  of  p52  NFkB2  (data  not  shown)  or 
the  NFkB  negative  regulator  IkBo.  PI3K-dependent  signal¬ 
ing  can  also  activate  NFkB  and  Akt  activation,  a  primary 
downstream  target  of  PI3K,  has  been  implicated  in 
antiestrogen  resistance.  However,  we  found  no  differences 
the  levels  of  phospho-Akt,  indicating  that  this  pathway  also 
is  unlikely  to  account  for  the  increased  NFkB  activity. 

In  contrast,  MCF7/LCC9  cells  express  ~  2-fold  higher 
levels  of  NEMO/IKKy.  NEMO  binds  to  IKKp  and  controls 
the  formation  of  the  IKK  complex  (37);  this  is  required  for 
the  activation  of  NFkB  in  response  to  external  stimuli  such 
as  tumor  necrosis  factor  a  (31,  38).  Up-regulation  of  NEMO 


Table  1.  Faslodex  and  parthenolide  synergistically  enhance 
apoptosis  in  MCF7/LCC9  cells 


Cell  line/drug 

%  Apoptosis  ±  SE 

P 

LCC1  vehicle 

4.22  ±  0.98 

— 

LCC1  ICI 

19.96  ±  4.43 

0.03* 

LCC9  vehicle 

3.20  ±  1.96 

— 

LCC9  ICI 

4.41  ±  0.90 

0.61* 

LCC9  parthenolide 

9.95  ±  1.21 

0.04* 

LCC9  ICI  +  parthenolide 

18.34  ±  1.45 

0.003*,  0.001  + ,  0.01  * 

NOTE:  RI  =  2.28  for  combination  of  1CI  and  parthenolide. 
^Versus  vehicle. 
tVersus  ICI. 
iVersus  parthenolide. 
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Figure  5.  Combined  treatment  with  Faslodex  and  parthenolide  has  no 
effect  on  the  stability  of  IkBo  expression.  LCC9  cells  were  treated  with 
100  nmol/L  Faslodex,  600  nmol/L  parthenolide,  Faslodex  +  parthenolide, 
or  ethanol  vehicle  in  CCS-IMEM  for  24  h  before  cell  lysis.  Lysates  (20  pg) 
were  separated  by  SDS-PAGE,  transferred  to  nitrocellulose  membrane, 
and  immunoblotted.  f. i-actin ,  loading  control. 

in  MCF7/LCC9  cells  would  enhance  the  kinase  activity  of 
IKK  and  likely  adds  to  the  elevated  levels  of  p65  to  further 
increase  basal  NFkB  activation. 

Constitutive  NFkB  activity  is  known  to  arise  as  breast 
cancer  cells  progress  to  an  estrogen-independent  (26,  27) 
and  antiestrogen-resistant  state  (10).  However,  this  is  the 
first  report  implicating  NEMO/IKKy  in  these  events. 
Regulatory  control  of  NEMO  is  complex,  involving 
sequential  small  ubiquitin-like  modifier  and  ubiquitin 
modification  occurring  in  both  the  cytoplasm  and  nucleus 
(39).  Whether  the  hormonal  regulation  of  NEMO  is  altered 
in  the  MCF7/LCC9  cells  has  not  been  determined,  but  is 
currently  being  pursued  to  clarify  further  the  mechanism 
by  which  NFkB  activity  is  elevated  in  breast  cancer  cells 
with  acquired  antiestrogen  resistance. 

The  NFkB  inhibitor  parthenolide  strongly  represses  the 
proliferation  of  MCF7/LCC9  cells  (100  nmol/L,  ineffective; 
ICso=  600  nmol/L)  and  restores  their  sensitivity  to  Faslodex. 
For  example,  whereas  treatment  with  100  nmol/L  Faslodex 
alone  is  ineffective,  50%  growth  inhibition  is  achieved  in  the 
presence  of  only  100  nmol/L  parthenolide.  This  interaction 
between  Faslodex  and  parthenolide,  which  generates  an 
estimated  RI  =  1.82,  is  synergistic  (29).  The  restoration  of 
Faslodex  sensitivity  by  parthenolide  is  a  significant  finding 
and  directly  supports  our  hypothesis  that  the  up-regulated 
NFkB  activity  present  in  MCF7/LCC9  cells  is  a  major 
contributor  to  the  antiestrogen  resistance  phenotype. 

We  have  explored  several  mechanisms  through  which 
parthenolide  may  restore  Faslodex  sensitivity  in  antiestro¬ 
gen-resistant  cells.  For  example,  NFkB  inhibition  could 
rescue  the  ability  of  Faslodex  to  block  ER-dependent 
transcriptional  activity.  Expression  of  several  well-charac¬ 
terized  estrogen-regulated  genes,  including  progesterone 
receptor,  pS2,  and  cnthepsin  D  (10),  is  increased  in  MCF7/ 
LCC9  cells.  This  may  reflect  the  7-fold  higher  basal  levels  of 
ER-dependent  transcription  in  vehicle-treated  MCF7/LCC9 
cells,  relative  to  that  seen  in  MCF7/LCC1  cells.  Faslodex 
completely  inhibits  transcription  of  an  ERE-luciferase 
construct  in  MCF7/LCC1  cells  (antiestrogen  sensitive)  but 
has  no  effect  on  ERE-mediated  transcription  in  MCF7/ 
LCC9  cells  (antiestrogen  resistant).  Whereas  exogenous 
expression  of  p65/RelA  NFkB  can  repress  ER  transcrip¬ 
tional  activity  in  vitro  (40),  treatment  with  parthenolide 
does  not  restore  the  ability  of  Faslodex  to  affect  ER 


transcriptional  function  in  MCF7/LCC9  cells.  Thus,  the 
elevated  basal  ER-dependent  transcriptional  activity  in 
these  cells  is  likely  not  due  to  NFkB  up-regulation. 
Parthenolide  may  rescue  traditional  Faslodex-mediated 
effects  in  these  cells  downstream  of  ER  activity  and/or  in 
nonclassic  ER  pathways. 

Faslodex  induces  G0-G[  arrest  in  MCF7/LCC1  cells 
treated  with  100  nmol/L  Faslodex  (data  not  shown),  but 
neither  Faslodex  nor  parthenolide  has  any  effect  on  MCF7/ 
LCC9  cell  cycle  distribution.  Parthenolide-mediated  G2 
arrest  has  been  observed  at  10-fold  higher  concentrations 
than  were  used  in  this  study  (35);  therefore,  the  possibility 
that  MCF7/LCC9  cell  cycle  progression  is  inhibited  by 
much  higher  concentrations  of  parthenolide  cannot  be 
excluded.  Indeed,  several  studies  with  parthenolide  and 
other  sesquiterpene  lactones  used  1  to  10  gmol/L  or  greater 
concentrations  to  achieve  50%  inhibition  of  cell  growth, 
whereas  our  cells  required  up  to  10-fold  lower  concen¬ 
trations  to  achieve  the  same  results. 

Our  results  show  that  the  inhibition  of  cell  growth  by 
Faslodex  and  parthenolide  is  not  primarily  cytostatic  in 
nature.  In  marked  contrast,  a  combination  of  Faslodex  and 
parthenolide  synergistically  promotes  programmed  cell 
death  (Rl  =  2.28).  Importantly,  the  proportion  of  apoptotic 
cells  observed  in  the  presence  of  both  drugs  (19%)  is 
comparable  to  that  seen  when  MCF7/LCC1  cells  are 
treated  with  the  same  dose  of  Faslodex  (19.9%).  Partheno¬ 
lide  can  enhance  the  apoptotic  activities  of  taxanes  and 
retinoids  (19,  20),  and  we  now  show  that  it  can  also 
potentiate  the  death  of  antiestrogen-resistant  breast  cancer 
cells  by  restoring  their  sensitivity  to  Faslodex. 

DeGraffenried  et  al.  (41)  have  recently  reported  that 
NFkB  inhibition  by  parthenolide  increased  breast  cancer 
cell  sensitivity  to  tamoxifen.  However,  these  investigators 
used  MCF-7  cells  genetically  engineered  to  overexpress 
activated  Akt;  these  cells  exhibit  tamoxifen  resistance  and 
NFkB  activation  that  is  entirely  dependent  on  Akt- 
mediated  pathways.  In  our  cell  system,  which  was  derived 
by  selection  in  the  presence  of  Faslodex  rather  than  by 
genetic  engineering  of  the  cells,  NFkB  up-regulation  does 
not  correlate  with  enhanced  Akt  activity. 

Sesquiterpene  lactones  in  general,  and  parthenolide  in 
particular,  can  prevent  the  degradation  of  IkB,  block 
activation  of  IKK,  alkylate  cysteine-38  in  p65/RelA  to 
prevent  DNA  binding,  and  inhibit  inducible  nitric  oxide 
synthase  (36,  42,  43).  We  found  no  evidence  of  NBa 
stabilization  when  MCF7/LCC9  cells  were  treated  with 
parthenolide  either  in  the  absence  or  presence  of  Faslodex. 
Parthenolide  can  also  inhibit  p42/44  mitogen-activated 
protein  kinase  activity  (42)  but  we  observed  no  reduction  in 
the  levels  of  phospho-mitogen-activated  protein  kinase  in 
our  cells  upon  parthenolide  treatment  (data  not  shown). 
Interestingly,  Nakshatri  et  al.  (44)  have  recently  shown  that 
parthenolide  can  reverse  breast  cancer  cell  resistance  to 
tumor  necrosis  factor-related  apoptosis-inducing  ligand  by 
enhancing  the  activation  of  c-Jun  NH2-terminal  kinase. 
Whether  parthenolide-induced  c-Jun  NH2-terminal  kinase 
activity  plays  a  role  in  its  restoration  of  Faslodex  sensitivity 
in  MCF7/LCC9  cells  has  yet  to  be  determined. 
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40  Parthenolide  Restores  Antiestrogen-Induced  Apoptosis 


Our  studies  clearly  show  that  treatment  with  the  NFkB 
inhibitor  parthenolide  is  a  viable  approach  to  restoring 
Faslodex-induced  apoptosis  in  breast  cancer  cells  that  have 
acquired  resistance.  Several  preclinical  studies  have  shown 
that  parthenolide  also  is  effective  in  the  treatment  or  chemo- 
prevention  of  cancer  cell  growth  (45).  A  phase  I  study  of 
feverfew  in  cancer  patients  was  recently  completed  and 
reported  no  significant  toxicity  observed  at  the  doses  tested 
(21).  Other  direct  or  indirect  inhibitors  of  the  NFkB  path¬ 
way  also  show  promise  as  antiproliferative  agents  (16)  and 
include  some  nonsteroidal  anti-inflammatory  drugs,  anti¬ 
oxidants,  immunosuppressants,  proteasome  inhibitors,  and 
glucocorticoids.  Our  work  now  shows  that  inhibition  of 
NFkB  may  also  be  useful  for  the  treatment  of  ER-positive 
breast  cancers  that  have  acquired  resistance  to  antiestrogen 
therapy,  thus  restoring  the  activity  of  one  of  the  most  active 
and  least  toxic  modalities  available  in  the  treatment  of 
endocrine-dependent  breast  cancer. 
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